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ABSTRACT 

We provide a systematic measurement of the rest-frame UV continuum slope /3 over a wide range 
in redshift {z^2-6) and rest-frame UV luminosity (0.1L*^3 to 2Ll^^) to improve estimates of the 
SFR density at high redshift. We utilize the deep optical and infrared data (ACS/NICMOS) over the 
CDF-South and HDF-North GOODS fields, as weU as the UDF for our primary UBVi "dropout" 
Lyman Break Galaxy sample. We also use strong lensing clusters to identify a population of very low 
luminosity, high-redshift dropout galaxies. We correct the observed distributions for both selection 
biases and photometric scatter. We find that the f/F-continuum slope of the most luminous galaxies 
is substantially redder at z~2-4 than it is at z^5-6 (from ~ —2.4 at z ~ 6 to ~ —1.5 at z ~ 2). 
Lower luminosity galaxies are also found to be bluer than higher luminosity galaxies at z ^ 2.5 and 
z ~ 4. We do not find a large number of galaxies with /3's as red as —1 in our dropout selections at 
z ~ 4, and particularly at z > 5, even though such sources could be readily selected from our data 
(and also from Balmer Break Galaxy searches at z 4). This suggests that star- forming galaxies 
^ ^ 5 almost universally have very blue C/y-continuum slopes, and that there are not likely to 
be a substantial number of dust-obscured galaxies at z > 5 that are missed in "dropout" searches. 
Using the same relation between C/V^- continuum slope and dust extinction as has been found to be 
appropriate at both z ^ and z ^ 2, we estimate the average dust extinction of galaxies as a function 
of redshift and UV luminosity in a consistent way. As expected, we find that the estimated dust 
extinction increases substantially with cosmic time for the most UV luminous galaxies, but remains 
small (< 2x) at all times for lower luminosity galaxies. Because these same lower luminosity galaxies 
dominate the luminosity density in the UV continuum, the overall dust extinction correction remains 
modest at all redshifts and the evolution of this correction with redshift is only modest. We include 
the contribution from ULIRGs in our star formation rate density estimates at z ^ 2-6, but find that 
they contribute only ^20% of the total at z ~ 2.5 and <10% at z > 4. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Quantifying the star-formation rate (SFR) density at 
z > 3 is a challenging endeavor. While there are a wide 
variety of techniques to estimate this rate at z < 2 from 
light at different wavelengths (Condon 1992; Kennicutt 
1998; Madau et al. 1998; RanaUi et al. 2003; Yiiksel et 
al. 2008; Li 2008; Reddy et al. 2006; Yun et al. 2001; 
Reddy & Steidel 2004; Reddy et al. 2006; Erb et al. 
2006b), the situation becomes considerably more uncer- 
tain at redshift z > 3. This is because most of the tech- 
niques effective at low redshift rely upon light at wave- 
lengths that simply cannot be detected at high redshifts 
for all but the most exceptional sources. Determinations 
of the SFR density from x-ray or radio emission appear 
to be just possible a,t z ^ 2-4 by stacking a large number 
of sources in very deep integrations (e.g., Seibert et al. 
2002; Nandra et al. 2002; Reddy & Steidel 2004), but 
do not work at z > 5 (Lehmer et al. 2005; Carilli et al. 



2008). Similarly, while deep mid-IR data would seem to 
be quite effective in estimating the SFR at z ^ 2 (e.g., 
Reddy et al. 2006; Caputi et al. 2007), at z > 3 many 
of the key spectral features redshift to wavelengths not 
very amenable to study with current instrumentation. 

This leaves us with sparingly few techniques for esti- 
mating the SFR density at high redshift, most of which 
have only recently been developed. These techniques in- 
clude traditional approaches like a consideration of the 
UV light to estimate the SFR densities to newer ap- 
proaches that use gamma-ray burst (GRB) rate densities 
(e.g., Yiiksel et al. 2008; Li 2008) or possible detections 
of redshifted Ha (e.g., Chary et al. 2005) to make these 
estimates. 

Perhaps, the simplest and most practical of these tech- 
niques is to rely on the UV light emanating from the 
high-redshift galaxies themselves. Young stars emit a 
large fraction of their energy at these wavelengths, and 
it is very easy using current instrumentation to measure 
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TABLE 1 

Observational data used to select z ^ 2 — 6 U, B, V, and i dropouts and to measure 

THE C/y-CONTINUUM SLOPE. 
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^ 5cr point-source limiting magnitude in 0.6" diameter apertures. 

While NICMOS Jno observations are available over both Abell 1689 and Abell 1703, these 
clusters lack deep coverage at ~ 1.6/im and so we do not make use of the existing NICMOS data 
to measure f/V-continuum slopes /3 for z ~ 5 V-dropout or z ~ 6 i-dropout selections. 



this energy to very low SFR rates (i.e., ~0.1 M© yr ^ at 
z ~ 4 in the Hubble Ultra Deep Field [HUDF: Beckwith 
et al. 2006]). One of the most significant challenges in 
establishing the SFR from the UV data is estimating the 
absorption by dust. Since UV light is subject to substan- 
tial attenuation by dust, this entire issue is quite impor- 
tant. One of the most practical methods for determin- 
ing the dust attenuation is through the J/F-continuum 
slopes (3 {fx oc A''). Since these slopes have been shown 
to be well correlated with the dust extinction at 2; 
(Meurcr ct al. 1995: Meurer et al. 1997; Mcurcr et al. 
1999; Burgarella et al. 2005; Dale et al. 2007; Treyer et 
al. 2007; Salim et al. 2007), 2; ~ 1 (Laird et al. 2005), 
and z - 2 (Reddy ct al. 2006: Daddi ct al. 2004), we can 
use these slopes to estimate the dust obscuration at even 
higher redshift. 

Not surprisingly, there have already been a number 
of attempts to determine these slopes at high redshift 
{zr^2-6) from the available imaging data. Adelberger 
& Steidel (2000) and Meurer et al. (1999) examined the 
UV slopes for a sample of z ^ 3 U -dropouts and inferred 
modest (factor of ^^5) dust extinctions. Other groups 
(Lehnert & Bremer 2003; Ouchi ct al. 2004; Stanway 
et al. 2005; Yan et al. 2005; Bouwens et al. 2006) have 
looked at these slopes at higher redshifts (z ~ 4 — 6), 
again on a dropout selection. Broadly, it was found that 
higher redshift galaxies had bluer J/y-continuum slopes 
than lower redshift galaxies have, suggesting that the 
typical dust extinction at high redshift is significantly 
less than it is at lower redshift. There was also some evi- 
dence from previous analyses (e.g., Meurer et al. 1999 
at 2: ~ 2.5) that lower luminosity galaxies had bluer 
?71^-continuum slopes 6 and thus lower dust extinctions 
(though the emphasis in Meurer et al. (1999) was the 
correlation between (3 and the dust-corrected luminosity, 
not the observed luminosity). 

Regrettably, apart from these two trends (higher red- 
shift galaxies are bluer and lower luminosity galaxies arc 
bluer), it has been difficult to make quantitative state- 
ments about the distribution of C/y-continuum slopes or 
dust extinctions for galaxies at high redshift. In part, 
this can be attributed to the piecemeal way the UV- 
continuum slopes and dust extinction have been derived 



at high redshift, as different analyses have derived UV- 
continuum slopes based upon a wide variety of disparate 
high-rcdshift galaxy samples, using different techniques 
(not all of which are consistent). In part, this can also 
be attributed to the limited amounts of data available 
to each particular study. Fortunately, the situation has 
begun to change, and there is now a wide variety of HST 
data that are available to select high redshift galaxies 
and determine their slopes over a wide range in redshift 
and luminosity. Determining how these slopes also de- 
pend upon luminosity seems particularly relevant given 
the fact that these sources likely dominate the UV lumi- 
nosity density and perhaps star formation rate density at 
z > 2 (e.g., Bouwens et al. 2006: Sawicki & Thompson 
2006b; Bouwens ct al. 2007: Reddy et al. 2008: Yan & 
Windhorst 2004; Reddy & Steidel 2009). 

It seems clear that a comprehensive, systematic de- 
termination of the distribution of J7T^-continuum slopes 
is needed, both as a function of redshift and luminos- 
ity. The goal of the present paper is to provide just 
such an analysis, taking advantage of the considerable 
quantity of deep HST data to clarify our understand- 
ing of dust obscuration and, by consequence, the star 
formation rate density at high redshift. We select galax- 
ies for J/F-continuum slope measurements over a wide 
range in redshift using well-tested z 2.5 U, z ^ A B, 
z ^ 5 V, and z ^ 6 i dropout selections. We use the 
wide-area ACS Great Observatories Origins Deep Sur- 
vey (GOODS) fields (Giavafisco et al. 2004a) to quantify 
these slopes at bright magnitudes and extend our anal- 
yses to very low luminosities by examining very deep 
HST data sets like the HUDF, and also by examining 
dropouts behind massive galaxy clusters where the sub- 
stantial magnification factors permit us to consider very 
faint sources. 

The plan for this manuscript is as follows. In §2, 
we present the observational data set. In §3, we de- 
scribe our procedure for generating source catalogs from 
these data, selecting the dropout samples, and estimat- 
ing the C/V-continuum slopes while correcting for the 
effect of object selection and photometric scatter. Fi- 
nally, we discuss our results (§4), examine the likely im- 
plications of these results for the effective dust extinc- 
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tion and SFR density at z~2-6 (§5), and then include 
a summary (§6). Throughout this work, we will find 
it convenient to quote results in terms of the luminos- 
ity L*^^ Steidel et al. (1999) derived at z ~ 3, i.e., 
Miyoo.AS = —21.07, for consistency with previous work 
- though we note that the Steidel et al. (1999) LF results 
are now updated (MiyocAS = —20.97 ± 0.14: Reddy & 
Steidel 2009) but still consistent with the previous de- 
termination. Where necessary, we assume fio = 0.3, 
J7a = 0.7, Hq = 70km/s/Mpc. Although these param- 
eters are slightly different from those determined from 
the WMAP five-year results (Dunkley et al. 2009), they 
allow for convenient comparison with other recent results 
expressed in a similar manner. Unless otherwise stated, 
the J/y-continuum continuum slope (3 presented here all 
assume a 1600A-2300A basefine. We express all magni- 
tudes in the AB system (Oke & Gunn 1983). 

2. OBSERVATIONAL DATA 

The purpose of this study is to measure the UV- 
continuum slope over as wide a range in luminosity and 
redshift as is practical from current observations. Of 
course, establishing this distribution even over a limited 
range in redshift or luminosity can be a challenge, as 
it requires very deep imaging over a wide wavelength 
range, both to select galaxies through the Lyman break 
technique and to measure their UV-continuum slope red- 
ward of this break. It also requires sufficient search area 
to provide a sufficiently large samples of sources to effec- 
tively define the distribution of [/^-continuum slopes. 

At z ^ 2.5, we use the deep [/-band coverage avail- 
able over the WFPC2 Hubble Deep Field (HDF) North 
(Williams et al. 1996) and HDF South (Casertano et al. 
2000) to select z ~ 2.5 [/-dropouts. Those two deep 
WFPC2 fields are invaluable for defining the distribu- 
tion of C/y-continuum slopes at lower luminosities, but 
their small area limits their value for establishing the 
distribution at the higher luminosities. However, as we 
will see, the C/lZ-continuum slope distribution for bright 
z ~ 2.5 galaxies is similar to what Adelberger & Steidel 
(2000) and Reddy et al. (2008) derived previously. For 
our ^ 4 i?, ~ 5 V^, and z ~ 6 i-dropout samples, 
we utilize the ACS i?435V606*775-Z850 data available over 
the HUDF, HUDF05 (Oesch et al. 2007; Bouwens et al. 
2007), and GOODS fields. The wide-area GOODS fields 
are used for establishing the C/V-continuum slope dis- 
tribution for luminous galaxies while the deeper HUDF 
and HUDF05 fields are used for establishing the UV- 
continuum slope distribution for lower luminosity galax- 
ies. The reduction of these latter data was performed 
with "apsis" (Blakeslee et al. 2003) and is described in 
Bouwens et al. (2006, 2007). A summary of the proper- 
ties of these fields is given in Table [H 

Importantly, the HUDF, HUDF05, and GOODS fields 
also have ~ 70 arcmin^ of very deep near-IR coverage 
with NICMOS - reaching 5a fimits of > 26.5 (0.6"- 
diameter apertures). This coverage is important for mea- 
surements of the t/y-continuum slope /3 for galaxies in 
our z ~ 5 ^-dropout and z ~ 6 i-dropout selections. 
The NICMOS data are described in more detail in R.J. 
Bouwens et al. (2009, in prep) and are illustrated in Fig- 
ure [T] (see also Bouwens et al. 2008). These data were ob- 
tained as a result of the cumulative observations of many 
different HST programs (e.g., Dickinson 1999; Thompson 



et al. 1999; Thompson et al. 2005; Bouwens et al. 2005; 
Oesch et al. 2007; Riess et al. 2007; C. Conselice et al. 
2009, in prep). All the NICMOS data taken since 2002 
(when the NICMOS cryocooler went into operation) were 
reduced with our NICMOS pipeline "nicred.py" (Magee 
et al. 2007). 

Finally, we take advantage of the deep HST data over 
three massive galaxy clusters to identify a small number 
of lower luminosity star-forming galaxies at z 4 — 6 
behind these clusters. Because of lensing by the fore- 
ground cluster (magnifying faint z ~ 5-6 galaxies by fac- 
tors of > 4) , it is possible to select and to measure UV- 
continuum slopes for a small number of extremely low 
luminosity galaxies. The five galaxy clusters (MS1358, 
CL0024, AbeU 2218, Abell 1689, AbeU 1703) all have 
deep ACS -B4355475V555r625«7752;850 data and cover 55 
arcmin^ in total. These clusters also have 6 arcmin^ of 
deep JiioHien NICMOS data necessary to measure the 
J/y-continuum slope for z ^ 5 V dropout and z 6 
i-dropout galaxies behind these clusters. The ACS data 
available over these clusters were also reduced with "ap- 
sis" while the NICMOS data were reduced with "ni- 
cred.py". These reductions are described in more detail 
in Bouwens et al. (2009a) and Zheng et al. (2009). 

3. ANALYSIS 

In this section, we describe our procedure for determin- 
ing the distribution of [/V^-continuum slopes (3 from our 
observational data as a function of redshift and UV lu- 
minosity. We begin by discussing how we do photometry 
on our observational data and then select z 2-6 galaxies 
using various dropout criteria (§3.1-3.2). In §3.3, we cal- 
ibrate our optical/near-IR photometry and then derive 
t/y-continuum slopes from the measured colors (§3.4). 
By combining measurements of the J7U-continuum slope 
determinations from both very deep and wide-area data 
sets, we derive the t/U-continuum slope distribution over 
a wide range in luminosity (§3.5). In §3.6, we take ad- 
vantage of lensing from massive galaxy clusters to ex- 
amine the J/y-continuum slope f3 distribution for a set 
of lower luminosity galaxies. We correct the observed 
distributions of [/U-continuum slopes /? for the effect of 
object selection and photometric scatter (§3.7). Finally, 
we discuss the uncertainties that likely exist in our de- 
terminations (§3.8) and then give our final estimates of 
the C/F-continuum slope distributions in §3.9. 

3.1. Photometry 

We used SExtractor (Bertin & Arnouts 1996) to per- 
form object detection and photometry. Object detec- 
tion was performed using the square root of the ini" 
ages (Szalay et al. 1999: essentially a coaddition of the 
relevant images) constructed from all ACS images with 
coverage redward of the dropout band (e.g., V606*7752850 
in the case of a i?435-dropout selection) or the WFPC2 
images redward of the dropout band in the case of our 
U300 dropout selection (where the -8450^606-^814 bands are 
used). For clarity, the image is equal to 

^k{hix,y)/Nkf 

where /^(a;, y) is the intensity of image Ik at pixel {x, y), 
where is the RMS noise on that image, and where 
the index k runs over all the relevant images (i.e., those 
redward of the dropout band) . 



4 





HDF^N GOODS 



Fig. 1. — Illustration of the deep near-IR data available over the two GOODS fields. Deep near-IR data are essential for measurements 
of the UV-continuum slope of galaxies at ^ > 5 and valuable for a measurement of these slopes at 2 ~ 4. The red and dark red regions 
correspond to regions with deep and very deep NICMOS Jiio and -ffieo coverage, respectively (5(t depths > 26.5 and > 28 AB mag). The 
light orange regions have deep NICMOS Hieo-band coverage {5cr depths of > 26.5: H. Teplitz et al. 2009, in prep; C. Conselice et al. 2009, 
in prep) while the magenta regions have deep NICMOS Jno-band coverage (5a- depths of > 26.5: R.J. Bouwens et al. 2009, in prep). The 
blue and dark blue regions correspond to regions with deep and very deep optical ACS V606*775Z850 coverage, respectively (5cr depths of 
> 28 and > 29 AB mag). For more details on the NICMOS data over the two GOODS fields, see Bouwens et al. (2008) and R.J. Bouwens 
et al. (200?, in prep). 



Our procedure for measuring colors depends upon the 
PSF of the data we are using. For colors that include 
only the ACS or WFPC2 bands (where the PSF is much 
sharper than for NICMOS), we measure these colors in 
scalable apertures determined using a Kron (1980) factor 
of 1.2. When measuring colors that include the NICMOS 
near-IR bands, we first PSF match the ACS/WFPC2 
data to the NICMOS i?i6o-band. Then, we use the same 
apertures described above if the area in those apertures 
is > 0.38 arcsec^ (i.e., >0.7"-diameter apertures). This 
only applies for the largest sources in our selection. Oth- 
erwise (and in most cases), the colors are measured in 
fixed 0.7"-diametcr apertures. 

We correct the fluxes measured in the above apertures 
(used for color measurements) to total fluxes using a two 
step procedure. First, the individual flux measurements 
are corrected to a much larger scalable aperture (using 
a Kron factor of 2.5). These corrections are made on a 
source-by-source basis using the square root of the 
image (Szalay et al. 1999). Second, a correction is made 
to account for the light outside of these larger apertures 
and on the wings of the ACS Wide Field Camera (WFC) 
PSF (Sirianni et al. 2005). Typical corrections are ~ 
0.6 mag for the first step and ~ 0.1 — 0.2 mag for the 
second step. While there are modest uncertainties in 
these corrections (typically 0.2 mag), these corrections 
only affect the total magnitude measurements and do 
not affect the colors. 

3.2. Sample Selection 

We select star-forming galaxies over the redshift range 
z ~ 6 to z ~ 2 using the well-established dropout 
technique (e.g., Steidel et al. 1996; Bunker et al. 2003; 
Vanzella et al. 2006; Dow-Hygelund et al. 2007; Stanway 
et al. 2007; Vanzella et al. 2009). Our f/-dropout cri- 
terion is similar to those used by Bouwens et al. (2003) 
while our B,V, and i dropout criteria are almost identi- 



cal to those adopted by Bouwens et al. (2007) except that 
our T^-dropout criterion uses a sharper («775 — zgso) < 0.6 
cut to exclude dusty z ~ 1 interlopers. These criteria are 

(C/300 - ^450 > 1-1) A (S450 - ^814 < 1-5) 
A(t/300 - ^450 > 0.66(5450 - hu) + 1.1) 



for our [/-dropout sample 



{B, 



435 



^606 > 1.1) A (B435 



^606 > (^606 - 
A(T/606 



^850) + 1-1) 
- Z850 < 1-6) 



for 



our B-dropout sample, 

[{V&OQ - «775 > 0.9(^775 " ^sso)) V (^606 " «775 > 2))] A 

(Veoe - «775 > 1-2) A (1775 - zsso < 0.6) 
for our y-dropout sample, and 

(i775-;2850 > 1.3)A((l/606-*775 > 2 {S / N {Vqqq) < 2)) 

for our ?-dropout sample, where A and V represent the 
logical AND and OR symbols, respectively, and S/N 
represents the signal-to-noise. These criteria are illus- 
trated in two color diagrams in Figure [2] and are very 
similar to those considered by Giavalisco et al. (2004b) 
and Beckwith et al. (2006). 

To ensure that our selections did not contain any spu- 
rious sources, we required that sources in our U , B, V , 
and i dropout selections be >5.5(T detections in the /8i4, 
^850 J ^850 J and zg,5o bands, respectively. We inspected all 
sources in our dropout samples by eye to purge them of 
any obvious contaminants (diffraction spikes, irregular- 
ities in the background). Pointlike sources (SExtractor 
stellarity S/G >0.8) were also removed (< 2% of the 
relevant sources). 

3.3. Calibration of the NICMOS Magnitudes 
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Fig. 2. — Two-color diagrams showing the selection criteria (see §3.2) we use to find z ~ 2.5 (7, 2 ~ 4 B, 2: ~ 5 V, and z ^ 6 i dropouts in 
our search fields (upper left, upper right, lower left, and lower right panel, respectively). Sources that fall within the shaded regions would 
be included in our selections. The blue tracks show the expected colors of starbursts with different l/V-continuum slopes as a function 
of redshift, while the red lines show the colors of low-redshift interlopers (Coleman et al. 1980). The black points illustrate the range of 
colors that real sources in our search fields possess. To ensure that these plotted distributions of colors are as realistic as possible, we only 
include sources from our deepest fields (the HDFs for our (7-dropout selection and HUDF for our B, V, and i-dropout selection) and at 
relatively bright magnitudes (I < 26 for our [/-dropout selection, 2 < 26.5 for our B-dropout selection, and 2 < 27 for our V,i-dropout 
selections). The red arrows show the Calzetti et al. (2000) reddening vectors. These diagrams show that our selection criteria are effective 
in identifying high-redshift galaxies with C/V-continuum slopes (3 bluer than 0.5, but less effective in identifying galaxies redward of this 
limit (Figures |4] and [6] show the effective selection volume calculated to be available to galaxies with various (/V-continuum slopes). 
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TABLE 2 

Wavebands used to derive the UV continuum slope for 
individual galaxies in high-redshift u, b, v, and i dropout 

SAMPLES.*''' 



Dropout 


Mean 


Filters used for 3 determination 


Sample 


Redshift 


Lower 


Upper 


(7-dropout 


2.5 


^606 (1730A) 


/814 (2330A) 


B-dropout 


3.8 


2775 (I6IOA) 


2850 (1940A) 


V-dropout 


5.0 


2850 (I55OA) 


(JllO + i?i60)/2 (2250A) 


V-dropout 


5.0 


2850 (I55OA) 


Hieo (2660A) 


i- dropout 


5.9 


Jiio (1570A) 


ifieo (2290A) 



* Sec §3.4 for discussion on the filter choices. 

^ Throughout this work, our ,fl measurements are based upon the 
relationship fxocX^. Since real spectra will not have perfect power 
law SEDs (and therefore depend upon the rest-frame wavelengths 
used to estimate /3), small corrections are made to the estimated /3's 
based upon a fiducial starburst SED (e^*/^ star formation history 
with t = 70 Myr, r = 10 Myr, [Z/Zq] = -0.7: from Papovich 
et al. 2001). See §3.4 and Appendix A for details. The filters that 
define the wavelength range of the /? measurements for our different 
dropout samples arc presented here. 

To accurately estimate the ?7y-continuum slopes /3 for 
galaxies in our sample, it is absolutely essential that the 
colors that we measure are accurate. For example, for 
a source with UV colors measured across the baseline 
I6OOA to 23OOA, a 0.1 mag error in this color translates 
into an error of 0.3 in the estimated slope /3. 

Measuring colors to this level of accuracy is challeng- 
ing for a number of reasons, particularly in the near-IR. 
First and foremost, the NICMOS detector (essential for 
measuring the fluxes of faint sources in the near-IR) has 
been shown to suffer from significant non-linearity in its 
count rate. While there are now procedures (e.g., de 
Jong et al. 2006) to correct for these non-linearities, as 
utilized here, it is unlikely that the corrections are per- 
fect in the magnitude range we are considering. This is 
because there are comparatively few deep fields available 
where the near-IR fluxes can be measured and calibrated 
against other data. Secondly, the colors we are measur- 
ing are derived from both optical (ACS) and near-IR 
(NICMOS) data where the PSF is very different. While 
such differences can be effectively addressed by carefully 
matching the PSFs of the two data sets before measuring 
fluxes (as we do here), PSF matching can be challenging 
to perform perfectly and so it is easy to introduce small 
systematic errors at this stage. 

As a result of these issues, we took great care in de- 
riving correction factors that we could apply to the mea- 
sured NICMOS near-IR fluxes. For simplicity, we as- 
sumed that these correction factors did not depend on 
the luminosity of the sources we were considering - but 
rather were simple offsets we could add to the measured 
magnitudes. We established these correction factors by 
examining the photometry of a sample (^50) of relatively 
faint (22 < zs5o ab < 26) point-like stars (SExtractor 
stellarity >0.8) over the HUDF/GOODS fields, and then 
compared their observed colors with that expected from 
the Pickles et al. (1998) atlas of stellar spectra. Given 
the very limited range in SED shape found in real stars, 
these comparisons should allow us to make reasonable 
corrections to the J and iJ-band photometry. We deter- 
mined the average differences between the measured and 



best-fit near-IR fluxes (both in the Jug and Hiqq bands) 
and found corrections of —0.13 mag and —0.07 mag for 
the Jiio and -ffieo band, respectively. 

To ensure that these corrections were as accurate as 
possible, we also identified a sample of ~ 60 highest 
S/N S dropouts from our search fields (GOODS and 
HUDF) and then fit their observed SEDs with Bruzual 
& Chariot (2003) r models. Since the 1.1/j, J and 1.6/j, 
H band fluxes of these galaxies probe the light blueward 
of the age sensitive break at 3600A (probing 2300A and 
3300A rest- frame, respectively), it should be reasonable 
to use the available i (I6OOA rest-frame) and 2;-band 
(19OOA rest-frame) fluxes to establish a Z7V-continuum 
slope p and extrapolate to these redder wavelengths. For 
our fiducial model fit parameters, we assmned a er^l"^ 
star formation history (with f = 70 Myr, t = 10 Myr), 
a Salpeter IMF, and [Z/Zq] = -0.7, and let the dust 
extinction and redshift be free parameters (the fiducial 
model considered by Papovich et al. 2001 in their stellar 
population modelling of z 3 [/-dropouts). Compar- 
ing the predicted near-IR fiuxes (from the models) with 
the measured fiuxes, we computed corrections to the J 
and ff-band fiuxes. The corrections were — 0.08±0.02 
mag and 0.04±0.2 mag, respectively (the uncertainties 
were derived using other SED models [Appendix A] to 
determine these corrections). Because of the size of the 
uncertainties in the iJ-band corrections (from the z ^ 4 
i?-dropouts), we gave those results very low weight over- 
all for the il-band correction. 

Combining the constraints we derived from our sam- 
ple of stars and z ~ 4 LBGs - which were consistent 
within 0.05 mag we derived a —0.10 mag correction to 
the measured Jug-band magnitudes. For the NICMOS 
7?i6o-band, we derived a —0.05 mag correction from the 
above fits - using the SED fit results to stars as our 
primary constraint.^ The measured Jno and i^ieo-band 
fluxes were therefore somewhat too faint (before correc- 
tion). A brief discussion of the effect that uncertainties 
in the photometry/zero points could have on our derived 
C/y-continuum slopes /? is given in §3.8. 

3.4. Measuring the UV -eontinuum slope 

Our principal interest is to determine the rest-frame 
C/y-continuum slope /3 for the galaxies in our dropout 
samples. As noted in the introduction, the UV- 
continuum slope /3 specifles how the flux density of a 
galaxy varies with wavelength (i.e., f\ oc A^) in the UV- 
continuum (i.e., from -^ISOOA to ~3500A). 

We estimate the C/y-continuum slopes /3 from the 
broadband colors available for each galaxy in our sample. 
To make our estimates of these slopes as uniform as pos- 
sible, we selected filters at each redshift where the rest- 
frame UV wavelengths were as very close to 1650A and 
23OOA as possible. Of course, the actual effective wave- 
lengths of these filters will vary somewhat from sample 
to sample - depending on the passbands in which broad- 
band imaging is available for our different samples or the 
precise redshift of specific galaxies in our samples. Fortu- 
nately, the wavelength range 1650 A to 2300 A is a good 

^ Note that the corrections we derive here are very similar (< 
0.04 mag) to what we derive matching sources in the HUDF with 
both NICMOS and WFC3/IR data (e.g., Bouwens et al. 2009b). 



7 



match to the wavelength range used for studies of the 
dust extinction in galaxies at z ~ (e.g., Meurer et al. 
1995; Meurer et al. 1999; Burgarella et al. 2005). Ta- 
ble [2] summarizes the passbands we use to probe these 
wavelengths for each of our dropout samples. 

Some additional clarification may be helpful regarding 
specific filter choices. For our B dropout sample, for ex- 
ample, we use the 1775 — zgso colors for these estimates to 
take advantage of the much higher quality (deeper, wider 
area) ACS data available. For our F-dropout sample, we 
use ^850 — ('^iio + HiQo)/2 colors predominantly for these 
estimates, but we also make use of zgso — i?i60 colors 
for these estimates to take advantage of the large ar- 
eas within the GOODS fields which have only NICMOS 
-ffi6o-band coverage (primarily from the Conselice et al. 
2009, in prep, NICMOS program: see the light orange 
squares in Figure [1]). 

In principle, deriving the f/l^-continuum slopes f3 from 
the measured UV colors is straightforward. One simply 
takes power-law spectra f\ cx A^, calculates their colors 
at the mean redshift of each of our dropout samples, 
and then determines for which P the model spectrum 
reproduces the observed colors. In practice, however, 
since the spectrum of star-forming galaxies in the UV- 
continuum is not expected to be a perfect power law (as 
noted above), small differences in the wavelength range 
over which the [/1^-continuum slope /3 is measured will 
have an effect of the derived results. 

To correct for these passband and wavelength- 
dependent effects, we base our estimates of /3 on a small 
number of more realistic SEDs calculated from stellar 
population models (Bruzual & Chariot 2003). These 
model SEDs are calibrated to have C/V^-continuum slopes 
P of —2.2, —1.5, and —0.8 over the wavelength range 
and assume a e"*/"^ star formation history {t = 70 Myr, 
r = 10 Myr, [Z/Zq] — —0.7) with varying amounts of 
dust extinction calculated according to the Calzetti et al. 
(2000) law. The formulae we use to convert between the 
measured colors and the C/V^- continuum slope (3 are pre- 
sented in Appendix A. The formulae we use to do the con- 
version depend somewhat on the star formation history 
assumed and likely result in an additional uncertainty in 
the derived [/F-continuum slopes f3 oi ^ 0.1 — 0.2. 

3.5. Constructing the UV -continuum Slope Distribution 
from Multiple (Deep or Wide-Area) Dropout 
Selections 

In order to accurately quantify the distribution of UV- 
continuum slopes /3 over a wide range in luminosity, we 
alternatively used our deepest selections and our wide- 
area selections. We used the deepest selections to define 
the ?71^-continuum slope distribution at the lowest lumi- 
nosities to maximize the S /N on the derived colors (while 
minimizing the importance of selection effects). Mean- 
while, we used our widest-area selections to define this 
distribution at higher luminosities. This allows us to find 
a sufficient number of sources to map out the distribution 
of [/l/-continuum slopes (which is helpful for determin- 
ing the mean and la scatter). In Table |31 we provide 
a list of the fields we use to quantify the distribution of 
t/y-continuum slopes at a given luminosity and redshift. 

In Figure [31 we plot the ?7F-continuum slopes /3 ob- 
served for galaxies in our four dropout samples versus 
their UV luminosities. The mean f/l^-continuum slope 



TABLE 3 

Dropout samples used to measure the distribution of 

C/y-CONTINUUM slopes /3 AS A FUNCTION OF REDSHIFT AND UV 
LUMINOSITY. 



Dropout 




Luminosity 


# of 


Sample 


Field 


Range'' 


Sources 



fZ-dropout 
_B-dropout 

V-dropout 
j-dropout 



HDF-North 
HDF-South 
GOODS-North 
GOODS-South 
HUDF 
Cluster 
GOODS-North 
GOODS-South 
HUDF 
HUDF05 
Cluster 
GOODS-North 
GOODS-South 
HUDF 
HUDF05 
Cluster 



-22 < MuvAB < -18 
-22 < Muv,AB < -18 
-23 < MuvAB < -20.5 
-23 < MuvAB < -20.5 
-20.5 < Muv,AB < -17 
-23 < Muv,AB < 
-23 < Muv,AB < 
-23 < Muv.AB < 
-23 < Muv.AB < 
-23 < Muv.AB < 
-23 < Muv.AB < 



AB 



z 

< 



-23 < Muv 
-23 < Muv./ 

-23 < Muv.AB < 
-23 < Muv.AB < 
-23 < Muv.AB < 



-20.5 
-20.0 
-20.0 
-19.5 
-19.5 
-20.0 
-20.5 
-20.5 
-20.0 
-20.0 
-20.5 



97 

71 

218 

251 

470 

195 

22 

19 

10 

16 

12 

2 

2 

5 

3 

2 



^ Luminosity limits for 2 ~ 4 _B-dropout samples arc chosen to 
minimize the effects of photometry scatter and selection on the 
(yV-continuum slope distribution. Luminosity limits for 2 ~ 5 V- 
dropout and 2 ~ 6 i-dropout samples are set according to the depth 
of the near-IR NICMOS data. See §3.5. 

P and la for galaxies is also included on this figure 
(blue squares) as a function of luminosity for each of the 
dropout samples. We adopt finer 0.5 mag bins for our 
B dropout selections than the 1.0 mag bins we adopt 
for our U and V dropout selections or 2.0 mag bins we 
adopt for our i dropout selection. The width of the bins 
depends upon the number of sources present in each of 
our dropout selections. The UV luminosity as shown on 
this diagram is the geometric mean of the two UV lu- 
minosities used to establish the UV slope (see Table [2]) . 
Use of the geometric mean of the two UV luminosities 
is preferred here for examining the correlation of (3 with 
luminosity since it allows us to evaluate the correlation 
without introducing any artificial correlations. Had we, 
for example, elected to use the bluer bands for examin- 
ing this correlation, the ijy-continuum slope /3 would be 
bluer as a function of the bluer-band luminosity, simply 
as a result of our examining the relationship as a function 
of this quantity. A similar (but opposite) bias would be 
introduced, in examining the ijy-continuum slope (3 as 
a function of the redder band. 

A clear trend is seen towards bluer C/y-continuum 
slopes P at lower luminosities. This trend is similar 
to that found by Meurer et al. (1999) in their analy- 
ses of z ~ 2.5 {/-dropout galaxies in the HDF North and 
also noted by Overzier et al. (2008) in their analysis of 
z ~ 4 g-dropouts in the TN1338 field (Figure 4 from that 
work) . 

3.6. UV -continuum Slopes Derived for z 4 — 6 
Galaxies Gravitationally Lensed By Galaxy 
Clusters 

In §3.5, we considered galaxies from both wide-area 
and deep surveys to establish the distribution of UV- 
continuum slopes P over a wide range in luminosity. We 
can increase the number of sources in our lower lumi- 
nosity z ~ 4 — 6 samples, by considering gravitationally 
lensed sources behind high-redshift galaxy clusters. The 
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clusters under consideration include Abell 2218, MS1358, 
CL0024, Abell 1689, and Abell 1703, and all allow us to 
select very faint star-forming galaxies at z 4 — 6 from 
the available HST ACS data. 

These clusters substantially amplify the flux from dis- 
tant galaxies, making it possible to measure the UV- 
continuum slope of galaxies at lower luminosities than 
would otherwise be possible. Of course, in the case 
of S-dropout selections from very deep optical data 
like HUDF, we are able to reach to the same intrinsi- 
cally luminosities as we reach in the i?-dropout samples 
we compile from searches behind lensing clusters. We 
correct the luminosities of the dropout galaxies iden- 
tified behind these clusters by applying gravitational 
lensing models from the literature. For Abell 2218, 
MS1358, CL0024 (sometimes known as CL0024-f 1652), 
Abell 1689, and Abell 1703, we adopt the models con- 
structed by Elfasdottir et al. (2007), Franx et al. (1997), 
Jee et al. (2007), Limousin et al. (2007), and Limousin 
et al. (2008), respectively. Typical magnification factors 
for the z ^ 4 — 6 galaxies we find behind galaxy clus- 
ters are ~5-10. We do not include the C/y-continuum 
slope results for sources where the model magnification 
factors are substantially greater than 10 - due to sizeable 
model-dependent uncertainties in the actual magnifica- 
tion and hence their luminosities (the uncertainties in the 
magnification factor for sources with such high predicted 
magnifications can be very large, i.e., greater than fac- 
tors of ~3-4: see Appendix A of R.J. Bouwens et al. 2009, 
in prep). Even for sources where the predicted magnifi- 
cation factors are smaller, the model magnification fac- 
tors (and hence the luminosities corrected for lensing) are 
likely uncertain at the factor of ~2 level (R.J. Bouwens 
et al. 2009, in prep). 

We include the C/y-continuum slope determinations of 
these galaxies with those estimated from our field sam- 
ples on Figure [3] (black points on z ~ 4 panel). As a 
result of a substantial magnification from gravitational 
lensing (typical factors of '^5-10), the intrinsic luminosi- 
ties of galaxies found behind clusters are much lower on 
average than those found in our field samples. In general, 
we observe good agreement between the distribution of 
i7F-continuum slopes /3 derived from our field samples 
and those inferred from our selections around clusters. 

3.7. Selection Biases and Measurement Uncertainties 

The distribution of C/y-continuum slopes that we de- 
rive is quite clearly affected by the manner in which 
sources are selected. Dropout criteria include galaxies 
with bluer L/y-continuum slopes more efficiently (i.e., 
over a larger range in redshift) than they do for galaxies 
with other colors. It is much easier to identify a sharp 
break in the SED of a blue galaxy than it is for a galaxy 
that is somewhat redder. For galaxies with red enough 
colors (i.e., f/^-continuum slopes /3 larger than 0.5), it 
is essentially impossible to robustly select high-redshift 
galaxy using the dropout criteria, and in fact the only 
sources that would satisfy the selection criteria would do 
so because of photometric scatter. This is illustrated in 
Figures [2] and HI 

To control for this effect, we construct models of the 
C/y-continuum slope /3 distribution, use these models to 
add artificial galaxies to real data, select sources from 
these data, and measure their [/l/-continuum slopes in 



the same way as from the real data. The goal is to 
construct a model that when "observed" reproduces the 
distribution of f/V^-continuum slopes /3 measured from 
the data. For the sizes and morphologies of the model 
galaxies used in the simulations, we start with the pixel- 
by-pixel profiles of the z ~ 4 i?-dropout sample from 
the HUDF (Bouwens et al. 2007) and scale their sizes as 
(1 -I- z)~^ (for fixed luminosity) to match the observed 
size-redshift trends (Bouwens et al. 2006; see also Fergu- 
son et al. 2004 and Bouwens et al. 2004). We use our 
well-tested "cloning" software (Bouwens et al. 1998a, b; 
Bouwens et al. 2003; Bouwens et al. 2006, 2007) to per- 
form these simulations. For the model UV LF at z ~ 2-6, 
we adopt the Schechter parameterizations determined by 
Reddy & Steidel (2009) at z ^ 2.5 and by Bouwens et 
al. (2007) at z ~ 4 - 6. 

In performing these simulations, we account for the 
modest correlation between the [/^-continuum slopes P 
of galaxies and their surface brightnesses. We determined 
the approximation correlation by examining 192 lumi- 
nous (~l-2 LIjy) galaxies in our GOODS _B-dropout se- 
lection and comparing their observed /3's with their sizes 
(half-light radii). Any correlation is potentially impor- 
tant, since it could substantially lower the efficiency with 
which we can select galaxies with very red f/y-continuum 
slopes P (e.g., Figure |4]). We find that sources with 
/9 > — 1 are ~15% larger than galaxies with — 2 < /3 < — 1 
and ^40% larger than galaxies (3 < —2 - though there is 
significant object-to-object scatter (the correlation coef- 
ficient between size (half-light radius) and /3 is just ~0.3). 
Fortunately, this correlation seems to only have a modest 
effect on these selection volumes, decreasing it by only 
~ 10% for the red galaxies and increasing it by only 
^10% for the blue galaxies. 

We experimented with a range of model UV- 
continuum slope P distributions to determine the effect of 
object selection and photometric scatter on the observed 
distribution of f/y-continuum slopes (3. These exper- 
iments were perfomed as a function of magnitude and 
the input color distribution (with mean t/y-continuum 
slopes ranging from —1.5 and —2.2 and the input color 
distribution taken to be Gaussian). In general, we found 
that the recovered distribution of /7V-continuum slopes 
P (after selection) is bluer than the input distributions 
(by A/3 ~ 0.1) at lower luminosities for all of our dropout 
samples. As expected, we found that noise in the ob- 
servations broadened the distribution of UV-contirmvLm 
slopes P somewhat over that present in the input distri- 
bution. Figure [4] illustrates how observational selection 
and noise modifies the input distribution of slopes for a 
i?-dropout selection over the HUDF. 

We used the simulations described above to quantify 
changes in this distribution as simple shifts in the mean 
C/U-continuum slope and la scatter. We used the esti- 
mated shifts to correct the observed distribution of UV- 
continuum slopes P for these effects. This corrected dis- 
tribution of L/U-continuum slopes is plotted in Figure [3] 
as open red squares, with the la scatter shown with the 
red error bars. 

To verify that the corrections we apply in this sec- 
tion are accurate, we compared the mean ?7y-continuum 
slope P we estimate near the faint-end (z85o_^b'^26.5- 
27) of our shallower GOODS _B-dropout selections (af- 
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Fig. 3. — Determinations of the l/V-continuum slope /3 (1600A to 2300A) versus rest-frame UV luminosity for the present selection of 
z ~ 2.5 C/, 2 ~ 4 B, 2 ~ 5 V, and z ~ 6 i dropouts. The UV luminosity here is the geometric mean of the UV luminosities measured in the 
bands used to establish the UV slope (see Table [2]l. On the right axes are the equivalent E{B — V) extinctions (Calzetti et al. 2000 law) 
to a given f3 if the base spectrum is a young star-forming galaxy (~100 Myr of constant star formation) - which is only AE{B — V) = 0.02 
bluer than our fiducial t model (i.e., the one referred to in footnote b of Table [2]l. The small points represent the C/y-continuum slopes 
measured for individual galaxies in our samples while the open blue squares and vertical bars show the mean t/y-continuum slope /3 and 
1(T scatter (shown in a darker blue in the 2 ~ 4 panel for clarity). Blue upper limits are shown for a few 2 ~ 5 and 2 ~ 6 sources with /3's 
bluer than —3.5. The green points on the 2 5 panel show the (yV-continuum slope determinations derived from our V-dropout samples 
where only NICMOS Hiao-hand data are available (i.e., with no deep NICMOS Jno-band coverage). The black points on the 2 ~ 4 panel 
show the {/V-continuum slope determinations derived from our B-dropout samples behind lensing clusters (where we have corrected their 
luminosities for model lensing magnification factors: see §3.6). It is encouraging that the {/V-continuum slopes we derive from both of the 
latter V and B dropout selections agree with the determinations derived from our primary selections. The open red squares and vertical 
bars show the mean C/V-continuum slope and Icr error in this slope, after correcting for the effect of object selection and photometric 
scatter (§3.7). See Table |4]for a tabulation of these slopes. The red lines show the best fit relationship between UV continuum slope and 
the observed magnitude of a galaxy in the UV (see §3.9 and §4.5). It is also clear that the (JV-continuum slope /3 of star-forming galaxies 
is much redder at higher luminosities than it is at lower luminosities (particularly at 2 ~ 2.5 and z ~ 4). 



ter correction) with those derived from our much higher 
S /N HUDF i3-dropout selections in the same magnitude 
range. We find that they are in excellent agreement, i.e., 
-1.81 ±0.04 at -19.2 AB mag for the GOODS fields vs. 
-1.73 ± 0.07 for the HUDF. This suggests that the cor- 
rections we apply in this section are reasonably accurate 
(A/3 < 0.1). 

3.8. Uncertainties and Model Dependencies 

Before presenting the distribution of {/^-continuum 
slopes /3 derived for each of the dropout selections, it is 
worthwhile to ask ourselves how our determinations may 
be affected by specific assumptions we make. We have 
already discussed the effect that errors in our photome- 
try (and zeropoints) would have on the ?7T^-continuum 
slope determinations (§3.3). A ^ 0.05 mag error in the 
derived colors would result in a 0.15 change in the de- 
rived f/T^-continuum slope /3. In an attempt to minimize 



such errors, we have exercised great care in obtaining a 
consistent set of colors across the optical and near-IR 
passbands for which we have data. 

We would expect a similar error in the [/V^-continuum 
slope P from the fiducial SEDs we use to convert from 
the observed colors to continuum slope /3. The un- 
certainties on the derived slope from this conversion are 
<0.1 for our z ~ 2.5, z 5, and z ^ 6 and <0.2 for our 
z ~ 4 sample (see Appendix A). 

The C/y-continuum slopes we derive for our selections 
also show some dependence upon the model redshift dis- 
tributions. As discussed in Appendix A, small differences 
(Az ^ 0.1) between the mean redshifts of dropout galax- 
ies in our models and that present in reality have a small 
effect on the derived C7T^- continuum slopes /3. Errors 
of Az ^ 0.1 in the mean redshifts of these selections 
shift the derived /3's by 0.02, 0.02, 0.01, and 0.05 for our 
z ^ 2.5 U, z ^ 4 B, z ^ 5 V, and z ^ 6 i dropout 
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Fig. 4. — {left) Comparison between the (/V-continuum slope /3 distribution input into the simulations (black histogram) and that directly 
recovered from the simulations {red histogram: which included the effect of object selection and measurement error). This panel is for 
our faintest, lowest S/N 2: ~ 4 _B-dropout selection over the HUDF (zgso-band magnitudes of 28-29). Both the input (and recovered) 
distributions are normalized to have the same area. On the top axes are the E{B — V) extinctions (Calzetti et al. 2000 law) that are 
equivalent to a given /3 if the base spectrum is a young star-forming galaxy (~100 Myr of constant star formation). The median UV- 
continuum slope /B input into the simulations and that recovered are shown with the black and red arrows, respectively. The recovered 
distribution of (71/-continuum slopes is broader and peaked to somewhat bluer values than the input distribution. For our higher S/N 
selections, the differences between the input and recovered distributions are much smaller (the offsets between the red and blue squares 
in Figure |3] show the size of the corrections). For our faintest selections (as shown here) the differences can be larger. It is nonetheless 
reassuring that when these corrections are applied there is agreement between the mean C/V-continuum slopes estimated in our shallower 
selections (where the corrections are larger) and in our deeper selections (where the corrections are smaller). For example, the mean 
C/V-continuum slope derived at 2g5o_AS~26.5-27 for our GOODS B-dropout selections (i.e., —1.81 ± 0.04 at —19.2 AB mag) agrees very 
well with that derived from our much higher S/N HUDF selections for the GOODS fields (i.e., -1.73 ± 0.07 for the HUDF). {right) The 
selection volumes (arbitrary units) available to galaxies for our B-dropout selection over the HUDF versus f/V-continuum slope /3. Shown 
are the relative selection volumes for galaxies with zgso ~24.7-26.7 (red line), zg^g ~26.7-27.7 (green line), and zgso ~27. 7-29.0 (blue line). 
We emphasize that the same normalization for the selection volume is used for all the three magnitude selections, demonstrating that it 
is much easier to select galaxies with red (7y-continuum slopes at bright magnitudes than it is at faint magnitudes. The reason bright 
2 ~ 4 galaxies are easier to select is that one needs to confirm that a source has a large Lyman break at >la, and this is easier to do at 
bright magnitudes. Selection biases and measurement errors are less for brighter samples (where the S/N is higher) and greater for fainter 
samples. The two panels here demonstrate how important it is for us to correct for selection biases and measurement uncertainties as we 
do in §3.7. 



selections, respectively. These are rauch smaller uncer- 
tainties than the A/3 ~ 0.15 errors that would result from 
small systematics in the photometry (or uncertainties in 
the conversion from observed colors to ?7y-continuum 
slopes). 

Similarly, the presence or absence of Lya emission (at 
1216A) in the spectra of sources in our selections is not 
expected to have a large effect on the measured UV- 
continuum slopes /3 themselves. This is because the mea- 
surements are made in passbands which are only sensitive 
to light redward of 1216A. The only exception to this is 
for our z ^ 6 i-dropout selection, but even there for typ- 
ical Lya EWs (~50A rest-frame: Dow-Hygelund et al. 
2007; Vanzeha et al. 2006; Stanway et al. 2007; Vanzella 
et al. 2009), the Jno — -ffieo colors would only change 
by ~0.07 - which would make the derived f3's ~ 0.17 
steeper. 

Nonetheless, the amount of flux in Lya has an effect on 
the redshift distribution of the sources selected with our 
dropout criteria. For galaxies with Lya EWs towards 



the upper end of the observed range (^50A rest-frame: 
Shapley et al. 2003; Dow-Hygelund et al. 2007; Vanzella 
et al. 2006; Stanway et al. 2007; Vanzella et al. 2009), 
the mean redshift of our dropout selections increases by 
Az ^ 0.3. This changes the derived /3's by 0.06, 0.06, 
0.03, and 0.15, respectively, for our four selections (Ap- 
pendix A). However, since only a small fraction (< 50%) 
of the star-forming galaxies at z ~ 2-6 appear to show 
Lya emission at these levels (Shapley et al. 2003; Dow- 
Hygelund et al. 2007; Vanzella et al. 2006; Stanway et 
al. 2007), the effect of the quoted uncertainties on (3 is 
likely much smaller than this (i.e., <0.075). Again, this 
is much smaller than the errors we would expect to result 
from small systematics in the photometry (or conversions 
to ?7T^-continuum slopes (3). 

Lastly, one might ask whether interstellar absorption 
features may have an effect on the t/y-continuum slopes 
P estimated from the broadband photometry. Fortu- 
nately, these absorption lines appear not to have a big 
effect on the measured slopes (i.e., A/3 < 0.1) as demon- 
strated by Meurer et al. (1999: see §3.4 and Figure 3 
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Fig. 5. — (top) Determinations of tlie mean l/V-continuum slope 
/3 for galaxies with a UV luminosity of L*_^ (Muy ~ —21) as 
a function of redshift (solid red circles). The error bars in /? are 
> 0.15 and include possible systematic errors in both the mea- 
sured colors and corrections for selection effects. The error bars in 
f3 (with the wider horizontal bars) are for the random errors. The 
random errors are generally much smaller than those including the 
estimated systematic errors (shown with the narrower horizontal 
bars). The right axes give the equivalent E{B — V) extinctions 
(Calzetti et al. 2000 law) for a given /3 if the input spectrum is a 
young star-forming galaxy (~100 Myr of constant star formation). 
Published determinations of the mean {/V-continuum slope /3 at 
2 ~ 3 (black star: Adelberger & Steidel 2000), at 2 ~ 4 (blue 
open pentagon: Ouchi et al. 2004), at 2 ^ 6 {solid cyan square: 
Stanway et al. 2005 and solid magenta triangle: Bouwens et al. 
2006), and at 2 ~ 4 — 6 {green open squares: Hathi et al. 2008) are 
also plotted. i*_3 galaxies are found to have bluer colors at 2 ~ 5-6 
than they do at 2 ~ 2-4. (bottom) Determinations of the mean UV- 
continuum slope /3 for galaxies with a UV luminosity of 0.15L*_3 
{Mjjy ^ —19) as a function of redshift (solid red squares). For 
comparison, our determinations of the mean (yV-continuum slope 
/3 at L*_g (from top panel) are included in this panel as solid red 
circles (and offset somewhat from the lower luminosity determina- 
tions in redshift for clarity). For each of our dropout samples, the 
lower luminosity galaxies are found to be bluer. 

from that work) given the wide wavelength range of the 
broadband filters used to estimate the slopes and the fact 
that most of these absorption lines occur at <1700A. 

3.9. Results 

In Table 21 we present the mean [/T^-continuum slopes 
and la scatter observed for our four dropout samples af- 
ter correction for selection effects and photometric scat- 
ter. This distribution is the same as that presented in 
Figure [3]). We assume a minimum systematic error in the 
mean f/y-continuum slope /3 of 0.15 as a result of pos- 
sible (~ 0.05 mag) systematics in the measured colors 
(see §3.3) and small model dependencies in the conver- 
sion from the observed colors to J/y-continuum slope /3 
(see Appendix A). 

There are clear trends that seem to be present in the 
distributions of J/lZ-continuum slopes /3 as a function of 
redshift and UV luminosity. The first trend is a corre- 
lation between the mean [/^-continuum slope and UV 
luminosity, in the sense that galaxies become bluer at 
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TABLE 4 

The mean f/y-CONTINUUM SLOPE /3 AND 1(7 
SCATTER OF GALAXIES, AS A FUNCTION OF UV 
LUMINOSITY.'*''' 



< Muv,AB > 



C/V-continuum slope (3 
Mean'^ la Scatter'* 



-21.73 
-20.73 
-19.73 
-18.73 

-22.22 
-21.72 
-21.22 
-20.72 
-20.22 
-19.72 
-19.22 
-18.72 
-18.22 
-17.72 
-17.22 

-21.90 
-20.90 
-19.90 
-18.90 

-20.76 
-18.76 



~ 2.5 (/-dropouts 
-1.18±0.17±0.15 
-1.58±0.10±0.15 
-1.54±0.06±0.15 
-1.88±0.05±0.15 

2 ~ 4 B-dropouts 
-1.32±0.12±0.15 
-1.33±0.07±0.15 
-1.44±0.04±0.15 
-1.51±0.03±0.15 
-1.58±0.10±0.15 
-1.68±0.09±0.15 
-1.73±0.07±0.15 
-1.93±0.06±0.15 
-1.98±0.05±0.15 
-2.03±0.04±0.15 
-1.91±0.05±0.15 

2 ~ 5 V-dropouts 
-1.81±0.22±0.15 
-1.76±0.10±0.15 
-2.17±0.08±0.15 
-2.64±0.16±0.15 

2; ~ 6 i-dropouts 
-2.16±0.15±0.15 
-2.32±0.19±0.15 



0.32 
0.27 
0.34 
0.36 

0.51 
0.31 
0.44 
0.40 
0.34 
0.35 
0.55 
0.31 
0.54« 
0.29'= 
0.26" 

0.45 
0.64 
0.66'= 
1.02" 

0.41 

0.91" 



" The mean C/V-continuum slopes /3 (I6OOA to 
23OOA) presented here are also shown on Figure|3] 
as the red squares. 

^ The slopes presented here have been corrected 
for selection effects and measurement errors (see 
§3.7). The tabulated UV luminosities are the ge- 
ometric mean of the UV luminosity measured in 
the two bands used to establish the UV slope (see 
Table [2J. 

Both random and systematic errors are quoted 
(presented first and second, respectively). In 
§3.8, we estimate the likely size of the system- 
atic errors. 

The la scatter presented here has been cor- 
rected for photometric scatter using the simula- 
tions described in §3.7 (see also Figure |4]l and 
therefore should reflect the intrinsic la scatter in 
the (/V-continuum slope (3 distribution. Typical 
uncertainties are ~0. 05-0. 10. 

" Because the observed scatter in the UV- 
continuum slopes for the faintest sources is dom- 
inated by the photometric errors, it is very dif- 
ficult to estimate the intrinsic la scatter in the 
f3 distribution, and therefore the uncertainties on 
our estimates of the intrinsic scatter in the /3 dis- 
tribution are large, i.e., > 0.2 — 0.4. 

lower UV luminosities. This trend is particularly signif- 
icant for our z ~ 2.5 t/-dropout and z 4 B-dropout 
selections, and fitting the mean /3 vs. UV luminosity 
to a line, we find P = (-0.20 ± OM){Muv.ab + 21) - 
(1.40±0.07±0.15) for our z ^ 2.5 [/-dropout sample and 
P = (-0.15±0.01)(M[/v,AB + 21)-(1.48±0.02±0.15)for 
our z ^ 4 B-dropout sample {the fit is shown on Figure\^ 
with the red lines)? The A/3 errors of ±0.15 given in the 
equations above are our estimates of the systematic er- 

^ While we find a strong correlation between the (/V-continuum 
slope /3 and the observed UV magnitude at 2 ~ 2.5 over the range 
Vsoe.AS ^ 23-27, it appears that the C/y-continuum slope /3 shows 
a weaker dependence on magnitude at brighter magnitudes, i.e., 
Veoe.AB ~ 23-25.5. Adelberger & Steidel (2000) and Reddy et al. 
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rors (see discussion in previous paragraph). We find no 
significant trend in tfie widtfi of the J/y-continuum slope 
P distribution as a function of luminosity. 

We also find a clear correlation between the mean UV- 
continuum slope /3 and redshift, in the sense that higher 
redshift galaxies are bluer. This is most evident in Fig- 
ure [5] where the mean UV-contiiumm. slope P is shown 
at a luminosity of L*^3 (i.e., M^y ~ -21) and 0.15L*^3 
(i.e., M^jy ^ —19) for each of our dropout samples. At 
z~2-4, the mean [/^-continuum slope is ~ —1.5 while 
at z > 5, it is < —1.8. No change is evident in the 
width of the J/y-continuum slope distribution - though 
this width is difficult to quantify for our highest redshift 
samples due to the small number of sources in these high 
redshift samples and the large photometric errors of the 
sources that are available. 

4. THE iyy-CONTINUUM SLOPE DISTRIBUTION: DIRECT 
IMPLICATIONS 

In §3, the C/y-continuum slope (3 distribution was de- 
rived as a function of both redshift and luminosity. Be- 
ing able to establish this distribution as a function of 
these two quantities is an important empirical result and 
has a number of noteworthy implications, which we will 
detail in this section. However, before describing these 
implications, it is prudent to compare the present UV- 
continuum slope (3 determinations with those in the lit- 
erature to put them in context (see §4.1). 

One of the most significant implications of these results 
is for the completeness of dropout selections at z ~ 2 — 6, 
which we discuss in §4.2. Being able to establish these 
distributions is also important for a determination of the 
z ~ 2 — 6 LFs (§4.3). In §4.4, we consider the question 
of how variations in the J/y-continuum slope /3 likely 
arise, and we argue that changes in the dust extinc- 
tion likely dominate the observed variations. Finally, in 
§4.5, we attempt to connect the sequence we find in the 
[7T^-continuum slope versus luminosity to similar trends 
found at lower redshift (and over other wavelength base- 
lines) . 

4.1. Comparison to previous determinations of the 
UV -continuum slope 

In §3, we derived the distribution of [/^-continuum 
slopes (3 over a wide range in redshift and UV luminos- 
ity using a very systematic approach while taking ad- 
vantage of a wide-variety of both deep and wide-area 
HST data. We presented evidence that the mean UV- 
continuum slope (3 is bluer at z ~ 5-6 than it is at z ~ 2-4 
and that this slope is also bluer at lower UV luminosities 
at z 2-4. 

Previously, there had been a variety of different at- 
tempts to measure these slopes at specific redshifts or 
luminosities (typically ~ ^2=3), e.g., Steidel et al. (1999), 
Meurer et al. (1999), Adelberger & Steidel (2000), Ouchi 
et al. (2004), Stanway et al. (2005), Bouwens et al. 
(2006), and Hathi et al. (2008). The top panel of Figure[5] 
provides a summary of many of these previous measure- 
ments. Most of the early high-redshift work focused on 

(2008) find no correlation between these quantities over this magni- 
tude range, and a t-test applied to our z ~ 2.5 f7-dropout sample 
shows a correlation at only 75% confidence over this magnitude 
range. 



z ^ 2-3 and was based upon U dropout selections from 
the HDF and large ground-based LBG searches (e.g., 
Steidel et al. 1999; Adelberger & Steidel 2000; Meurer 
et al. 1999). In those papers, the f/lZ-continuum slope p 
was found to have a mean value of ^ —1.4 and la dis- 
persion of ^^0.5-0.6 at ^ L^jy. No significant correlation 
of f/V^-continuum slope (3 with UV magnitude was found 
to -20 AB mag (Adelberger & Steidel 2000; Reddy et al. 
2008), though there would appear to be some evidence 
in the /? distributions presented by Meurer et al. (1999: 
e.g.. Figure 5 from that paper) that the (3 distribution be- 
comes a little bluer at lower UV luminosities (i.e., —18.5 
AB mag).^ 

Ouchi et al. (2004) extended these studies to higher 
redshift by presenting a determination of the UV- 
continuum slope distribution at 2 ~ 4 based on a large 
selection of B dropout galaxies from the Subaru Deep 
Field and Subaru XMM/Newton Deep Field. Ouchi et 
al. (2004) found that the C/lZ-continuum slope distribu- 
tion at z ~ 4 was consistent with that determined at 
z ^ 3 and that there was only a marginal trend (not 
statistically significant) towards bluer slopes at lower lu- 
minosities. Papovich ct al. (2004), by contrast, working 
with a selection of z ~ 4 i? dropouts from the GOODS 
fields, found that galaxies at z ~ 4 were slightly bluer in 
their UV colors than at z ^ 2.5 in the Hubble Deep 
Field North and South. At somewhat fainter magni- 
tudes, Beckwith et al. (2006) remarked that the UV 
colors of z ~ 4 _B dropouts in the HUDF were very 
blue in general (with /3's of ~— 2) and hence suggested 
very little dust extinction. The somewhat bluer UV- 
continuum slopes found by Beckwith et al. (2006) than 
by e.g. Ouchi et al. (2004) is not particularly surprising 
given the substantial differences in the mean luminosities 
of the two samples. Incidentally, the trend towards bluer 
J/y-continuum slopes at lower luminosities reported on 
here is evident in Figure 18 of Beckwith et al. (2006) 
although it was not specifically noted. 

At even higher redshifts, Lehnert & Bremer (2003) 
found that z ^ b V dropouts had t/y-continuum slopes 
P very close to —2 while later Stanway et al. (2005) and 
Bouwens et al. (2006) found [/^-continuum slopes of 
—2.2 ± 0.2 and —2.0 ± 0.3, respectively, at z '~ 6 from a 
selection of i-dropouts in the HUDF. Hathi et al. (2008) 
used a small sample of galaxies at z 4 — 6 from the 
HUDF and larger sample of [/-dropouts at z ~ 3 in an 
attempt to quantify the change in the f/y-continuum 
slope P as a. function of redshift. Hathi et al. (2008) 
found that the [7V^- continuum slope P showed only a 
slight change from z '--^ 3 (where P ~ —1.5) to z^5-6 
(where P 1.8). 

Broadly, the picture that has emerged from these stud- 
ies is that star-forming galaxies become bluer towards 
higher redshifts and to lower luminosities, but it has been 
difficult, given possible systematics between the different 
studies, to quantify the size of the changes. We confirm 
this overall picture, finding that the mean [/T^-continuum 
slope is ^ 1.0±0.4 bluer at z 6 than it is at z '--^ 2.5 and 
~ 0.5 ± 0.1 bluer at -18 AB mag (- 0.1^*^3) than it is 

Unfortunately, Meurer et al. (1999) do not provide a lot of dis- 
cussion on the possible correlation of /3 with observed UV luminos- 
ity (despite the existence of a likely trend) and instead emphasizes 
the correlation of /3 with dust-corrected UV luminosity. 
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at —21 AB mag (~ L^^^). Because of the larger samples 
and use of a consistent approach in deriving these slopes 
at all redshifts and luminosities, the differential evolution 
we measure is more robust than in previous studies. 

Looking more specifically at the [/V-continuum slope f3 
measurements we have derived at various redshifts and 
comparing those measurements with those obtained in 
previous studies, we find very good agreement in gen- 
eral, in almost all cases within the quoted errors. The 
only significant exception to this is the Hathi et al. (2008) 
measurements of the [/y-continuum slope /3 at z ~ 6 
where a mean /3 of —1.73 ± 0.35 was found. This ap- 
pears to result from the relation Hathi et al. (2008) use 
to compute the [/^-continuum slope /3 at z ~ 6 (i.e., 
(3 = 2.56(J — H) — 2). This relation does not account 
for the fact that the NICMOS Jno band extends down 
to 8OOOA and therefore z ~ 6 i-dropouts partially drop 
out in the Jhq band (i.e., are fainter in the Jno band 
and hence have redder Jno — Hjmo colors). Comparing 
P = 2.56(J - H)-2 with Eg. 1^ from Appendix A, we 
can readily see why the mean f/V^-continuum slopes we 
derive are much steeper (by ~0.3). We note that an ad- 
ditional —0.1 mag shift in /? relative to previous 2; ~ 6 
measurements of /3 (e.g., Bouwens et al. 2006; Hathi et al. 
2008) comes from small offsets we make to the Jno-band 
photometry (§3.3). 

4.2. Do LBG selections miss a substantial population of 
red star-forming galaxies at z > A? 

The Lyman Break galaxy selection technique identi- 
fies galaxies at high redshift through simple color crite- 
ria. This technique is well established to be an efficient 
and robust method for identifying star-forming galaxies 
over a wide range in redshift z~2-6 (Steidel et al. 1996; 
Williams et al. 1996; Bunker et al. 2003; Vanzella et al. 
2006; Dow-Hygelund et al. 2007; Vanzella et al. 2009). 
This technique provides a very complete census of UV 
light at high redshift - simply by virtue of the selection 
wavelength itself. However, it is less efficient for prob- 
ing the total stellar mass or even the total SFR at high 
redshift (e.g., van Dokkum et al. 2006). This is because 
galaxies with the highest stellar masses or SFRs are of- 
ten either old or dust obscured - making them fainter in 
the UV and thus more difficult to select with the LBG 
technique. 

Our interest here is in examining the systemic com- 
pleteness of LBG selections. We want to determine 
whether there is a set of galaxies at high redshift that 
we miss altogether by virtue of the LBG selection tech- 
nique. Note that this is a very different question from 
determining whether there is a class of galaxy at high 
redshift that we select less efficiently because they are 
faint in the UV (e.g., because of dust or age). In general, 
we would expect sources to miss our LBG selections if 
one of the two LBG color criteria failed to hold, i.e., if 
(1) the sources did not show a strong Lyman break or (2) 
the sources were too red in their i/y-continuum slope to 
satisfy the LBG selection. We would expect criterion (1) 
to always hold for a sufficiently high redshift source, as a 
result of line blanketing by the Lya forest. However, we 
might expect criterion (2) to fail if the ?7T^-continuum 
slopes f3 of the sources were sufficiently red. 

We can attempt to address this question by looking 
at how many galaxies have t/y-continuum slopes /3 that 



lie close to the selection limits. If our samples contain 
a large number of such galaxies, it would suggest that 
our LBG selections suffer from significant incompleteness 
near these limits. Figure [2] shows the selection criteria for 
each of LBG selections and Figure [S] presents the effec- 
tive search volume (§3.7) for galaxies {dashed black lines 
on the left-hand panels) as a function of ?7y-continuum 
slope (3 (see also Figure |4]). It is quite clear from these 
figures that our selection criteria is effective in identify- 
ing star-forming galaxies at 2: 2-6 with /3's bluer than 
^ 0.5, and it is striking how much redder this limit is 
than the [/V^-continuum slopes derived for galaxies in 
our dropout samples at z~2-6 {solid blue histograms in 
Figure [6]). The situation is particularly conspicuous for 
galaxies with t/y-continuum slopes /3 redder than —1 at 
z > 5. While our simulations (§3.7) suggest that galaxies 
with these colors should show up in our selections if they 
existed {dashed lines in Figure (HI), essentially none are 
found. This suggests that star-forming galaxies with red 
C/y-continuum slopes (3 are extremely rare. 

Complementary Balmer Break Selections: Another 
way we can investigate the issue of possible systemic in- 
completeness in LBG selections is through other selec- 
tion techniques. One such technique identifies galaxies 
based upon their rest-frame optical light and searches 
for a prominent Balmer Break at ^3700A rest-frame. 
Such selections are not surprisingly called Balmer Break 
Galaxy (BBG) selections and typically identify the old- 
est and most massive galaxies. Since such galaxies tend 
to be more chemically evolved, they are frequently more 
dusty. Considering such selections therefore permit us 
to evaluate the extent to which our LBG selections may 
miss redder and more dust-obscured starbursts. Bram- 
mer & van Dokkum (2007) performed such a selection 
at 2<2:<3 and 3<z<4.5 (see right-hand panels on Fig- 
ure [6]) based upon the Faint Infrared Extragalactic Sur- 
vey (FIRES) data (Labbe et al. 2003; Forster Schreiber 
et al. 2006) and derived J/y-continuum slopes for sources 
in their selections from SED fits to the photometry. They 
found that ~67% of the galaxies in their z ^ 2-3 selection 
had C/y-continuum slopes bluer than 0.5, but almost all 
(> 90%) of the galaxies in their z ~3-4.5 sample had 
such blue slopes. This is highly encouraging for 2 ~ 4 
LBG selections given the independent nature of Balmer 
Break selections - suggesting that z ~ 4 LBG selections 
may be largely complete. 

By contrast, the fact that ~33% of the galaxies in the 
Brammer & van Dokkum (2007) z 2-3 BBG sample 
have /3's redder than 0.5 suggests that completeness could 
be somewhat of a concern for LBG selections at z ~ 2.5, 
and in fact it is well known that at z 2-3, there is a 
substantial population of dust-obscured galaxies under- 
going vigorous star formation (e.g., Hughes et al. 1998; 
Barger et al. 1999; Chapman et al. 2005; Labbe et al. 
2005; Papovich et al. 2006; Pope et al. 2006). 

Searches for Balmer Break galaxies at z > 5 have also 
been conducted (e.g., Dunlop et al. 2007; Wiklind et al. 
2008: see also Rodighiero et al. 2007 and Mancini et 
al. 2009). As with the Brammer & van Dokkum (2007) 
study, these searches would seem to be relevant to the 
present discussion we are having about the completeness 
of LBG selections for star-forming galaxies at high red- 
shift. Unfortunately, there is a wide diversity of search 
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Fig. 6. — (left) Distribution of l/V-continuum slopes {3 observed at 2:~2-6 for the present Lyman-Break Galaxy selections (solid blue 
histograms). We only include the observed (3 values for Lyman-Break galaxies >1 mag brighter than the selection limit to avoid being 
dominated by selection effects. On the top axes are the E{B — V) extinctions (Calzetti et al. 2000 law) that are equivalent to a given 
l3 if the base spectrum is a young star-forming galaxy 100 Myr of constant star formation). The selection volume (§3.7) available for 
galaxies with C/V-continuum slope /3 and magnitudes > 1 mag brighter than the selection limit is shown with the thick dashed lines in each 
panel (the normalization is arbitrary). While it is more difficult to select galaxies with redder (/V-continuum slopes /3 with a Lyman Break 
criterion, our simulations (§3.7) indicate that our dropout selections should be capable of identifying modest numbers of galaxies with /3's 
as red as 0.5. Indeed, such galaxies (with /3's redder than —1) are present in our 2~2-4 LBG selections albeit in small numbers. However, 
they are conspicuously absent in our 2^5 selections. This suggests that such sources are exceedingly rare at z > 5. Conversely, the modest 
number of galaxies with red t/V-continuum slopes /3's seen at 2 ~ 4 seem consistent with the small samples of 2 ~ 4 ULIRGs that have 
been found in the observations (e.g., Daddi et al. 2009). The lack of very blue /3 < —3 galaxies is also consistent with expectations from 
models and suggest that the /3 distributions derived here are reliable, (right) Distribution of (yV-continuum slopes observed for Lyman 
Break Galaxy selections at 2 ~ 2-4 (solid blue histogram) and for Balmer-Break Galaxy selections at 2 ~ 2-4 (red histogram: Brammer 
& van Dokkum 2007). The E{B — V) extinctions given on the top axis are as in the left panel (but are not expected to be particularly 
representative for the somewhat older BBGs). Both /3 distributions plotted are as observed, with no correction for selection or measurement 
effects. Balmer-Break Galaxy (BBG) selections provide us with an independent measure of the (JV-continuum slope distribution since they 
do not depend upon their C/y-continuum slope or the brightness of these sources in the (/V-continuum. The fact that the distribution 
of C/y-continuum slopes /3 is becoming substantially bluer for both selections (and in particular the Brammer & van Dokkum 2007 BBG 
selection at 2 ^ 4) suggests that Lyman-Break selections (sensitive to galaxies with /3's bluer than 0.5: see Figures [2] and should be 
largely complete at 2 > 5 (see §4.2). 
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results in this area that make it difiicuh to draw clear 
conclusions. Wiklind et al. (2008) report 11 plausible 
Balmer Break galaxies at 2: ~ 5-7 over the CDF South 
GOODS field, 7 of which are detected at 24)j, with MIPS, 
while Dunlop et al. (2007) find no Balmer Break galaxies 
over the same redshift range, in their analysis of the same 
field. While the MIPS detections for 7 of the 11 z~ 5- 
7 Wiklind et al. (2008) BBG candidates may be inter- 
preted as due to obscured AGN, we believe a much more 
likely explanation is due to PAH emission, as Chary et 
al. (2007b) argue for HUDF-JD2 (Mobasher et al. 2005). 
The other BBG candidates from Wiklind et al. (2008) 
may have redshifts of 2; 5-7, but may also have lower 
redshifts. 

Implications: Putting together the present LBG se- 
lection results at z 2-6 with those from Balmer Break 
selections over the same range and other results in the 
literature, we find clear evidence that the galaxy popula- 
tion at high redshift is increasingly blue as one moves out 
to high redshift, and therefore high-redshift LBG selec- 
tions seem likely to be increasingly complete (and con- 
sequently suffer much less from systematic incomplete- 
ness). These trends are very clear in LBG selections at 
z ~ 2-6, and since there is no reason to suppose that these 
trends come from various selection biases (not only do the 
sources have /3's very far from the selection limits, but 
the effect of observational selection are corrected for), it 
seems reasonable to take the observed evolution at face 
value. Supporting evidence comes from a similar evolu- 
tion seen in complementary Balmer Break selections (at 
least according to Brammer & van Dokkum 2007). Pre- 
viously, Bouwens et al. (2007: §4.1) have discussed this 
in the context of z ~ 4 i?-dropout selections. 

Obviously we would expect galaxies at very high red- 
shifts (z ~ 5-6) to have much bluer UV-contummn slopes 
than at z ~ 2 — 3, because of the much smaller time base- 
line over which to produce metals and dust, as well as 
the shorter dynamical time scales of galaxies at z > 5 
(and hence much reduced ages). We would also expect 
this result based upon the evolution of the UV LF. Since 
the characteristic luminosity of galaxies in the UV be- 
comes progressively smaller as we move to higher red- 
shifts (Dickinson et al. 2004; Shimasaku et al. 2005; 
Bouwens et al. 2006; Yoshida et al. 2006; Bouwens et 
al. 2007; Oesch et al. 2009), we would expect to find 
fewer and fewer galaxies at high redshift with very large 
SFRs. Since dust extinction is well correlated with the 
SFR (e.g., Wang & Heckman 1996; Hopkins et al. 2001; 
Martin et al. 2005; Reddy et al. 2006; Buat et al. 2007; 
Zheng et al. 2007), we would not expect to find many 
galaxies at z > 5 with substantial dust extinction. Fi- 
nally there is the evolution seen in the relationship be- 
tween SFR and dust extinction. From z to z ^ 0.7 to 
z ~ 2, it has been found that the effective dust extinction 
for a given star formation rate decreases quite strongly as 
we move out to higher redshifts (Reddy et al. 2006; Buat 
et al. 2007; Burgarella et al. 2007). Each of these con- 
siderations suggest that galaxies at high-redshift should 
be almost uniformally very blue and LBG selections very 
complete. 

Nonetheless, we know there is a substantial population 
of luminous, dust obscured galaxies at z ^ 2 — 3 (e.g., 
Hughes et al. 1998; Barger et al. 1999; Chapman et al. 
2005; Labbe et al. 2005; Papovich et al. 2006; Pope et al. 



2006) and extending out to redshifts as high as z ~ 4.5 
(e.g., Capak et al. 2008; Daddi et al. 2009; Wang et al. 
2009; Dannerbauer et al. 2008). While some authors 
(e.g., Daddi et al. 2009) make the suggestion that these 
populations might be quite significant and contribute 
substantially to the overall SFR density at z > 4, we 
would argue that the role of these galaxies is likely much 
more modest in scope and they do not provide a large 
fraction of the SFR density, particularly at z > 4. We 
discuss and quantify this point in §6.2 (see also Tabled 
Table H Figure [H and Figure [HI). 

4.3. Relevance for LF Determinations 

As shown in Figure [2l Figure |4] and Figure [6] (and as 
discussed in the previous section), the effective selection 
volume for star-forming galaxies at z ~ 2-6 is very sen- 
sitive to its J/y-continuum slope (3 and UV luminosity 
and can vary by factors of > 2x (see also discussion in 
Sawicki & Thompson 2006a or Beckwith et al. 2006). 
Quantifying the distribution of C/y-continuum slopes for 
star-forming galaxies (as a function of both luminosity 
and redshift) is therefore an important first step in de- 
termining the UV LFs at high redshift. 

Given this impact, it has been somewhat surprising 
that different teams have made very different assump- 
tions about the distribution of J/y-continuum slopes in 
their determinations of the UV LFs. Some teams have 
adopted C/y-continuum slopes of f3 ^ —1.4 (Sawicki & 
Thompson 2006a; Yoshida et al. 2006), other teams have 
assumed UV continuum slopes of (3 ^ —2.0 (Beckwith et 
al. 2006), and yet other teams have adopted luminosity- 
dependent C/V^-continuum slopes (Bouwens et al. 2007). 
Such differences have only added to the large dispersion 
seen in LF determinations at high redshift (see Figures 
10, 11, and 13 from Bouwens et al. 2007). 

As a result of the present quantification of UV- 
continuum slopes (3 at high redshift (vs. redshift and 
luminosity), we now have a relatively uniform set of as- 
sumptions they can be used for quantifying the UV LFs 
at z ~ 2-6 and in extending such determinations to z > 7. 

4.4. Interpreting variations in the UV -continuum slope 

P 

Having used the available observations to derive the 
distribution of J/V-continuum slopes /3 for high-redshift 
galaxies as a function oi UV luminosity and a range in 
redshift (Table |4] and Figure [3|), we might ask ourselves 
what this teaches us about the physical properties of 
high-redshift galaxies. Since the [/T^-continuum slope 
13 is purely an observational measure of the shape of the 
spectrum of high-redshift galaxies and can be affected by 
a variety of different physical conditions (including dust, 
age, metallicity, etc.), we cannot use the above measure- 
ments to make any unambiguous observational inferences 
about the nature of high-redshift galaxies. 

Nevertheless, we will argue that the simplest way to ex- 
plain most of observed differences in [3 is through changes 
in the dust content of galaxies. While there are un- 
doubtably variations in other quantities such as the age, 
metallicity, or IMF of the stars that affect these slopes, 
we will argue that variations in the dust content of galax- 
ies likely have the largest effect on the observed UV- 
continuum slopes and we can use the observed variations 
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Fig. 7.— C/y-continuum slope /3 {I6OOA to 2300A rest-frame 
baseline) predicted for specific fractional changes in the mean age, 
metallicity, or dust extinction of high-redshift star-forming galax- 
ies from our fiducial stellar population model. For our fiducial 
model (where /3 = —1.47), we assume t = 70 Myr, r = 10 Myr, 
[ZIZq] = -0.7, E{B -V)= 0.15, and a Salpeter IMF (where 
the star formation history is parametrized as e~*/'^) from the Pa- 
povich et al. (2001) fits to 2 ~ 2.5 C/-dropouts from the WFPC2 
HDF North. In modifying our fiducial model to have younger 
ages, we make changes to both t and t. A factor of 2 (0.3 dex) 
change (dahsed line) in the mean dust content, age, or metallicity 
of the star-forming galaxy population at high redshift shifts the 
C/V-continuum slope /3 by ~ 0.35, ~ 0.1, and ^ 0.05, respectively. 
The effect of the stellar IMF on f3 is not shown here, but results 
in only modest changes in /3, i.e., A/3 < 0.1 for a ~ 0.5 shift in 
the slope of the IMF. It seems clear that changes in the mean dust 
content of galaxies at high-redshift has the biggest effect on the 
(yV-continuum slope /3 - though some change is also likely the re- 
sult of differences in the average age. Nonetheless, since stellar 
population modelling of 2; ~ 4 — 6 dropout galaxies (Stark et al. 
2009) suggest that age shows little dependence on redshift or lumi- 
nosity (at least in the median), it seems likely that age only plays 
a minor role (A/3 < 0.2) in driving the observed trends (§4.4). 

in these slopes to make inferences about changes in the 
effective dust extinction as a function of galaxy luminos- 
ity and redshift. 

We note that we would not expect the presence or lu- 
minosity of an AGN to have a big effect on these slopes, 
given that the observed incidence of such sources at z ~ 2- 
4 is just ^ 3% for moderately bright galaxies at z ~ 2 — 3 
(e.g., Nandra et al. 2002) and there is little evidence they 
are more frequent or important at higher redshift or lower 
luminosities (e.g., Lehmer et al. 2005; Ouchi et al. 2008). 

To ascertain which of the aforementioned factors (e.g., 
the overall dust content and properties, the age, metal- 
licity, or IMF of a stellar population) are likely to be the 
most important for interpreting the variations we see in 
the [/T^-continuum slopes (3, we consider the effect that 
changes in many of the above quantities would have on 
the C/y-continuum slope p. For simplicity, we model the 
star formation history of galaxies with a simple r model 
e"*/"^ with T = 10 Myr (e.g., as in Papovich et al. 2001) 
and t equal to the age. We also take the metallicity to be 
fixed and not evolve over this entire history. Finally, we 
implement the dust extinction by applying the Calzetti 
et al. (2000) law to the SED resulting from the stellar 



population models (Bruzual & Chariot 2003). Using the 
Papovich et al. (2001) modelling of z ~ 2.5 [/-dropouts 
in the HDF North as a guide, we adopt t — 70 Myr, 
Salpeter IMF, [Z/Zq] = -0.7, and E{B - V) = 0.15 as 
our fiducial parameters. 

We then change the age, metallicity, and dust content 
of this model by various factors and calculate the effect it 
would have on the [/V-continuum slope /3 predicted for 
a galaxy, as measured over the basehne I6OOA to 2300A. 
The results are shown in Figure [T] From this figure, it is 
clear that the largest changes in the C/y-continuum slope 
come from variations in the dust content of galaxies and 
the other parameters have a smaller effect. For factor of 
~2 changes in the age, metallicity, and dust, we estimate 
that the C/y-continuum slope (3 change by '^0.05, '^0.1, 
~0.35, respectively. 

Of course, we also see that age has a modest effect on 
the observed C/y-continuum slope f3. This is well docu- 
mented in the literature (e.g.. Bell 2002; Kong et al. 2004; 
Cortese et al. 2006; Panuzzo et al. 2007). Despite this 
sensitivity of f3 to age, no significant change are found in 
the median ages of star-forming galaxies from z ^ 6 to 
z ~ 4 (Stark et al. 2009), and little changes are found 
in the median ages as a function of luminosity (there 
is a hint that lower luminosity galaxies may be ~1.5x 
younger). Galaxies are found to have median ages of 
'^150 Myr. This suggests that the effect that systematic 
changes in the age of the stellar populations on the mean 
(3 is at most modest (< 0.2). Moreover, even if we ignore 
these observational results, it seems unlikely that the av- 
erage ages of star-forming galaxies would increase more 
rapidly than some multiple of the dynamical time in a 
galaxy. From theory (e.g.. Mo et al. 1998), the dynamical 
times scale as p^^l"^ (1 + z)^"^/^ - which is a factor of 
^ 3 (0.5 dex) from z 6 to z ^ 2.5. For such a change 
in age, the change in /? would be A/3 — 0.2, which is 
small compared to the observed differences - where A/3 
is 0.5-1.0. 

Metallicity has an even smaller effect (by a factor of 
~8) on the observed C/y-continuum slopes /3 than either 
dust or age do, so we would require very large variations 
in the metallicity to affect the C/y-continuum slope /3 in 
any sizeable way. However, since any substantial change 
in metallicity would almost certainly be accompanied by 
a similar change in the dust content (given the correlation 
between these two quantities), we would again be left 
with a situation where the changes in j3 resulting from 
metallicity would be completely overwhelmed by changes 
in /3 resulting from dust. 

Changing the IMF of galaxies only appears to have a 
modest effect on the C/y-continuum slope /3. For exam- 
ple, changing the slope of the IMF by ~0.5 only changes 
/3 by ;$0.1. Moreover, for a steep enough [/^-continuum 
slope and a top heavy enough IMF, there is very red neb- 
ular continuum emission (resulting from the ionizing ra- 
diation: Schaerer 2002; Venkatesan et al. 2003; Schaerer 
2003; Zackrisson et al. 2008; Schaerer & de Barros 2009) 
that more than offsets the very blue spectrum from the 
stars themselves (Figure 1 from Schaerer & Pello 2005). 

In summary, we expect that the observed variations in 
the mean [/^-continuum slope (3 to be largely the result 
of changes in the mean dust extinction. 
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4.5. Connection to trends found in star-forming 
galaxies z ~ — 3.' 

One of the most salient trends we found in the UV- 
continuum slopes (3 derived for z ~ 2.5 [/-dropout and 
z ^ 4 B-dropout galaxies was the presence of a strong 
correlation between the C/y-continuum slope (3 and UV 
luminosity (see also Meurer et al. 1999). As detailed in 
§3.9, the mean C/V-continuum slope j3 for z ~ 2.5 and 
z ^ A galaxies decreases (becomes bluer) by 0.20±0.04 
and 0.15±0.01, respectively, for each magnitude we reach 
fainter in luminosity. 

It seems reasonable to imagine that such trends might 
also be present in even lower redshift samples or in other 
colors for z^2-A samples. In fact, a strong correlation 
between [/T^-optical colors and rest-frame optical mag- 
nitude has been found for star-forming galaxies in the 
"blue" cloud (e.g., Papovich et al. 2001; Baldry et al. 
2004; Papovich et al. 2004; Wyder et al. 2007; Labbe 
et al. 2007), in the sense that more luminous galaxies 
are redder. This is very similar to the relationship we 
find between the [/T^-continuum slope (3 and UV lumi- 
nosity. Looking at the comparison more quantitatively, 
the slope of the color-magnitude relationship Labbe et al. 
(2007) measure, for example, is equivalent to dP/dMjjv 

of 0.27. This is similar to the -0.20 ± 0.04 slope 

we estimate at z ~ 2.5 for our [/-dropout selection (see 

Finally, it is worthwhile to remark on the physical ori- 
gin of this slope - which could plausibly be explained 
by changes in either the mean age or the dust content 
of galaxies as a function of luminosity. In the previous 
section, we argued that the simplest way of accommodat- 
ing the sizeable change (A/3 ^ 0.6) in the [/T^-continuum 
slope P with luminosity was through a change in the dust 
content of the galaxy population. Explaining the change 
in P with age would require very large changes (factor 
of ~10 changes in the mean galaxy age) - which seem 
contrary to the modest changes in age noted by Stark et 
al. (2009) as a function of luminosity. 

Labbe et al. (2007) also argue that the slope in the 
color-magnitude relationship is predominantly the result 
of a variation in the dust content. Labbe et al. (2007) 
draw this conclusion based upon an examination of the 
slope of the color-magnitude relationship vs. wavelength 
and through a detailed examination of the z ~ Nearby 
Galaxy Field Sample (Jansen et al. 2000). This sug- 
gests that the color-magnitude relationship we observe 
for star-forming galaxies may simply be another mani- 
festation of the well-known mass-metallicity relationship 
observed at z ~ — 3 (e.g., Tremonti et al. 2004; Erb et 
al. 2006a; Maiolino et al. 2008). 

5. INFERRED DUST EXTINCTION 

In §4.4-§4.5, we argued that the most likely inter- 
pretation of the systematic changes in the mean UV- 
continuum slope /? (as a function of redshift or luminos- 
ity) is a change in the dust content of galaxies. Given 
the sizeable (factor of ~5) dust corrections inferred at 
z ~ 2 — 3 for luminous galaxies (e.g., Reddy et al. 2006; 
Meurer et al. 1999; Erb et al. 2006b; Reddy & Stei- 
del 2004), any changes in the estimated dust corrections 
would have a substantial effect on the SFR densities in- 
ferred at earlier redshifts and hence its evolution across 



cosmic time. 

In this section, we use the measured [/F-continuum 
slope f3 distribution at z ~ 2-6 to estimate the dust cor- 
rections for [/IZ-bright galaxies at z~2-6. We begin by 
describing the formula we use to estimate dust extinction 
using the measured [/^-continuum slopes (§5.1). Then, 
we discuss the likely significance of the observed trends 
in [/y-continuum slope and dust extinction, versus lu- 
minosity (§5.2) and redshift (§5.3). In §5.4, we combine 
the derived trends in dust extinction with the observed 
UV LFs to calculate dust corrections, to various limiting 
luminosities. We then use the typical dust corrections 
for specific luminosity galaxies to estimate the bolomet- 
ric luminosities of specific luminosity galaxies in the UV 
(§5.5). Finally, in §5.6, we discuss the few well-cited cases 
of dusty high-redshift galaxies and explain why they do 
not alter the simple picture we lay out in this section. 

5.1. Inferred Dust Extinction 

To estimate the effective dust extinction for galaxies 
in our various samples and as a function of redshift, we 
rely upon the correlation found between dust extinction 
and the [/T^-continuum slope /3 at z ^ (Meurer et al. 
1999): 

Aieoo = 4.43 + 1.99/3' (1) 

where the dust extinction Aiqqq here is specified at rest- 
frame I6OOA and where /?' is the [/^-continuum slope 
measured from 1300 A to 2600 A. This relationship has 
been shown to be a reasonable predictor of the actual 
dust extinction in galaxies atz~0,z^l,Z'^2 for all 
but the youngest (Reddy et al. 2006; Siana et al. 2008; 
Siana et al. 2009) and most obscured starburst galaxies 
(e.g., Meurer et al. 1995; Meurer et al. 1999; Burgarella 
et al. 2005; Laird et al. 2005; Reddy et al. 2006), and 
therefore it is reasonable for us to use this equation to 
estimate extinction in each of our high-redshift dropout 
samples. This technique has already been employed in 
several previous studies estimating the SFR density at 
z~2-6 (e.g., Adelberger & Steidel 2000; Meurer et al. 
1999; Bouwens et al. 2006; Stark et al. 2007). 

Of course, there are many reasons for suspecting that 
Eq[l]may not work very well at early times, both because 
of the very different ages and metallicities of stellar pop- 
ulations and because the dust itself may have very dif- 
ferent characteristics. After all, the standard mechanism 
for forming dust in the envelopes of AGB stars should 
not work at z > 5 because the universe is not old enough 
at these redshifts to have produced AGB stars. Dust in 
2 > 5 galaxies has therefore been suggested to form by 
another mechanism, most frequently in SNe ejecta (e.g., 
Maiolino et al. 2004; Maiolino 2006). Eq.[T]is also known 
to fail at z 3 for a few lensed sources (e.g., MS1512- 
cB58 [Siana et al. 2008] or the Cosmic Eye [Siana et 
al. 2009]) and for very young (<100 Myr) star-forming 
galaxies at z ~ 2 (Reddy et al. 2006). 

To use Eq. [H we need to convert the [/V^-continuum 
slopes P we have derived (which have a baseline I6OOA to 
23OOA) to that appropriate for Eq. [T] (where /3 assumes 
a slightly more extended baseline 1300 A to 2600 A). We 
therefore explored the value of the [/^-continuum slope 
over these two different baselines for a wide variety of 
star formation histories and values of the dust extinc- 
tion. We found that they are very similar in general, with 
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TABLE 5 

The effective dust extinction (at ~1600A) 
estimated for the lbg population integrated 
down to various uv luminosities (see also 
Figure [sJ . 



Effective Extinction 



Sample 



> 



Using Meurer et al. (1999) Relationship' 



(7-dropouts 
S-dropouts 
V-dropouts 
i-dropouts 



6.0 
5.8 
2.7 
1. 



1.8+2.1c 
1.4-1.6 
+0.8+2.1 
0.7-1.5 
+0.7 + 1.0 
0.5-0.7 
+0.4+0.6 
-0.3-0.4 



Q O+0.8 + 1.4 
'-'■°-0. 6-1.0 

O+0.5+1.3 
■'■"-0.5-1.0 

1 (,+0.3+0.6 
^■"-0.2-0.4 

,+0.3+0.5 



1.3^ 



0.3-0.4 



2.8 
2 

1.4 
1.2 



a,b 

+ 0.5 + 1.0 
0.4-0.7 
+0.3+1.0 
-0.3-0.7 
+0.2+0.5 
0.2-0.4 
+0.2+0.4 
0.2-0.3 



The effective dust extinctions given here are the mul- 
tiplicative factors needed to correct the observed UV lu- 
minosity densities at ~ 1600 A to their intrinsic values, af- 
ter integrating to specific limiting luminosities (with the 
limits specified in the columns). These extinctions are 
estimated based upon the distribution of f/V-continuum 
slopes observed (Table |4jl and the correlation between 
(JV-continuum slope and dust obscuration observed at 
2^0-2 (Eq. [1] Meurer et al. 1999). Notice that the 
dust extinctions are much lower when integrated to very 
low luminosities. 

Both random and systematic errors are quoted (pre- 
sented first and second, respectively). 

This estimate is in good agreement with the estimates 
at z ~ 2 by Reddy et al. (2006) and Erb et al. (2006b) 
from various multiwavelength data (§5.4). 

an absolute difference generally A/3<0.2. While there is 
clearly latitude in the precise relationship we use to con- 
vert the measured C/y-continuum slopes (3 to this new 
wavelength, perhaps the least arbitrary set of SEDs to 
use to perform this conversion are those we use in Ap- 
pendix A to estimate /? from the observed colors. For 
those SEDs (which assume a e^*/^ star formation history, 
t = 10 Myr, T = 70 Myr, [Z/Zq] = -0.7, Calzetti et al. 
2000 extinction law, and a Salpeter IMF), we derive the 
following relationship /3' = -2.31 + l.ll(/3 + 2.3). While 
obviously there are some uncertainties in using this con- 
version, there are at least as many uncertainties in using 
Eq. [T] to estimate the dust obscuration in high redshift 
galaxies, and so we will ignore these uncertainties in the 
subsequent discussion. 

5.2. Dependence of Dust extinction on UV Luminosity 

Wc will first discuss the dependence of the UV- 
continuum slope and inferred dust extinction on the UV 
luminosity. We find that the ?7y-continuum slope /3 
shows a strong dependence on the UV luminosity at both 
z 2.5 and z ~ 4. This suggests that dust extinction is 
positively correlated with the observed UV luminosity, 
in the sense that more luminous galaxies are more dust 
obscured. 

At face value, this trend would appear to be con- 
trary to what Adelberger & Steidel (2000) and Reddy 
et al. (2008) found in their examination of z ~ 3 ?7- 
dropouts - where no correlation was found between the 
J/y-continuum slope /3 and the observed UV luminos- 
ity over the range 7?, = 22 {Muv,ab ~ —23.5) to 
TZ ^ 25.5 {MiJv,AB ~ —20). This is a somewhat brighter 
range than we probe with good statistics, suggesting that 
dust extinction (and the [/T^-continuum slope /3) may 
be a weaker function of the observed UV magnitude at 
Muv,AB ^ -20 (JZ < 25.5) than it is over the whole 
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Fig. 8. — The effective dust extinction (at ~1600A) estimated 
for LBGs vs. redshift integrated down to a luminosity of 0.3L*_3 
(red lines), 0.04L*_3 (green lines), and (blue lines). The precise 
redshifts at which our dust extinction estimates are presented are 
shifted slightly to minimize overlap between the different limiting 
luminosities. The error bars shown here are for the random errors. 
The dotted error bar included on the right here is the systematic er- 
ror that likely applies to each dust extinction estimate (based upon 
our conservative estimate that all /3 estimates in this paper are un- 
certain at the A/3 ~ ±0.15 level). This systematic error is unlikely 
to have any effect on the trends and differences seen here, but could 
change the overall scaling. The dust extinctions presented here are 
also given in Table |5] These extinctions are estimated based upon 
the distribution of C/V-continuum slopes observed (Table |4)l and 
the correlation between C/V-continuum slope and dust obscuration 
observed at 2 ~ - 2 (Eq. [T] Meurer et al. 1999). Notice that 
the dust extinctions are much lower when integrated to very low 
luminosities. 

range probed here. Indeed, using this same magnitude 
baseline (Vi306,AB < 25.5) and a t test, we only find the 
correlation between (3 and UV luminosity to be signifi- 
cant at 75% confidence from our own z ~ 2.5 [/-dropout 
sample. Thus, while the correlation would appear to be 
weaker for the more luminous sources, the existence of a 
correlation between (3 and observed UV magnitude over 
the entire magnitude range probed here is very signifi- 
cant (i.e., >4(T result). 

In order to understand the differences between our re- 
sults and those of Adelberger & Steidel (2000), we must 
appreciate that galaxies exhibit very different behaviors 
as a function of UV luminosity, depending upon whether 
we are dealing with galaxies at higher or lower luminosi- 
ties. For the lower luminosity regime, a large fraction of 
the sources should have very low star formation rates. In 
light of the well-known correlation between the star for- 
mation rate and dust extinction (e.g., Wang & Heckman 
1996; Hopkins et al. 2001; Martin et al. 2005; Reddy et 
al. 2006; Zheng et al. 2007), these galaxies would have 
low dust extinction (less than < 0.2 mag), and therefore 
the UV luminosity should be approximately proportional 
to the star formation rate. Consequently, we might ex- 
pect the dust extinction (and [/T^-continuum slope) to 
be approximately proportional to UV luminosity at low 
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luminosities. This is what we observe. 

In the higher luminosity regime, we might expect the 
SFR to be much larger in general. We would also expect 
the dust extinction in these galaxies to also be larger, 
given the correlation of dust extinction with the SFR. 
As such, any increase in the SFR of a galaxy would be 
largely offset by a similar increase in the dust extinction 
- resulting in a very weak correlation between the UV 
light output from a galaxy and its SFR (or dust extinc- 
tion). This is exactly what Adelberger & Steidel (2000) 
and Reddy et al. (2008) found in their sample of lumi- 
nous (> 0.3L*^3) 2^3 [/-dropout galaxies. Of course, 
we would also expect galaxies with very large extinctions 
to be present in galaxy samples at lower UV luminosi- 
ties, albeit as a small fraction. We can infer this latter 
fact from the volume density of such systems: > 10^^ 
Lq IR ultra- luminous galaxies (e.g., from the z ~ 2 IR 
LF of Caputi et al. 2007) have a much smaller volume 
density than UV faint sources (by factors of >50: see dis- 
cussion in §8 of Reddy & Steidel 2009). Consequently, 
we would expect the contribution from such IR ultra lu- 
minous galaxies to the overall SFR density in UV faint 
samples to be relatively small. 

5.3. Dependence of Dust Extinction on Redshift 

Very high-redshift {z > 5) galaxies are also much bluer 
in general than z 2-4 galaxies of the same luminosity. 
This suggests z > 5 galaxies are less dust obscured as 
well. This smaller dust obscuration is probably a conse- 
quence of an evolution in the relationship between dust 
extinction and the star formation rate, such that for 
a given star formation rate the observed dust obscura- 
tion increases as a function of cosmic time. This would 
cause higher redshift galaxies of a given UV luminosity 
to have bluer [/T^-continuum slopes. Such an evolution 
in extinction-SFR relationship is observed from z ~ 2 to 
z ~ (Reddy et al. 2006; Buat et al. 2007; Burgarella et 
al. 2007), and it is expected that this relationship con- 
tinues to evolve from z 6. This is presumably a result 
of the gradual build-up in metals in the universe with 
cosmic time. 

5.4. Average Dust Extinction for LBGs to Specific 
Luminosity Limits 

It is interesting to estimate what the observed distri- 
bution of C/y-continuum slopes (3 would imply for the 
dust extinction, if Eq. [T] holds. For this calculation, we 
assume that we can represent the distribution of UV- 
continuum slopes /? as a simple gaussian with mean and 
la scatter as derived earlier in our paper (Table ID see 
§3.7). We then fold this distribution through Eq. [1] to 
calculate the effective dust extinction integrated down 
to various limiting luminosities. We adopt the UV lumi- 
nosity function of Bouwens et al. (2007) at z ~ 4, z ^ 5, 
and z ~ 6 and that of Reddy & Steidel (2009) at z - 3 
for this calculation. 

The effective dust extinction is presented in Table [5] 
and Figure [S] To a limiting luminosity of 0.3 Ll^^, 
we infer a dust extinction of 6.0±2.5x (or 1.6 mag) at 
z ~ 2.5, very close to the 5.4±0.9x dust attention Meurer 
et al. (1999) estimated previously from the UV contin- 
uum slopes of HDF-North [/-dropouts. It is encourag- 
ing that this extinction estimate is very similar to the 
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Fig. 9. — Fraction of the total W luminosity density that arises 
from galaxies brighter than a given luminosity (solid line). This 
calculation assumes the faint-end slope a of —1.73 preferred in 
recent determinations of the LF at z ~ 2-6 (e.g., Yoshida et al. 2006; 
Bouwens et al. 2007; Oesch et al. 2007; Reddy k. Steidel 2009). The 
dashed lines show the fraction of the luminosity density brighter 
than specific luminosities for faint-end slopes ol of —1.68 and —1.78 
(which differ at the Icr level from the o = —1.73 value preferred 
at 2 ~ 4: Bouwens et al. 2007). Note that > 50% and > 80% of 
the luminosity density occurs in galaxies fainter than 10~'^ '^L* ~ 
0.05L* and 10~'''*-L* ~ 0.4L*, respectively, demonstrating how 
important it is to understand the dust correction down to very low 
luminosities. 



values estimated at z 2 using a wide variety of mul- 
tiwavelength data and SFR calibrators (e.g., Reddy & 
Steidel 2004; Reddy et al. 2006; Erb et al. 2006b). How- 
ever, these extinction estimates are a few times higher 
than those inferred by Carilli et al. (2008) at z ~ 3 by 
stacking the Very Large Array (VLA) radio observations 
of > 0.2L*^3 [/-dropouts in the COSMOS field (where 
a factor of 1.8±0.4 extinction was inferred). This may 
indicate that the Meurer et al. (1999) relationship sub- 
stantially overestimates the dust extinction at z > 2, or 
it may indicate that the radio luminosity for a given star 
formation rate is systematically lower at high redshift 
(Carilli et al. 2008). At this time, it is not obvious how 
this issue will be resolved, though it seems clear that 
dust extinction will be less of a concern as we reach back 
further in cosmic time (e.g., see §6.2). 

At higher redshifts and integrated down to lower lu- 
minosities, the effective dust extinction is much less, i.e., 
factors of <2-3 (< 1 mag). Qualitatively, we would ex- 
pect such a change due to the much bluer UV contin- 
uum slopes /? found at these luminosities and redshifts. 
Nonetheless it is still quite striking how much smaller the 
estimated dust extinction is when including light from 
the lowest luminosity galaxies. In detail, this is because 
lower luminosity galaxies dominate the total luminosity 
density in the UV continuum (see Figure (HI) and the av- 
erage dust extinction in these lower luminosity galaxies 
is quite low. Bouwens et al. (2007) estimated that half 
of the light in the [/T^-continuum originated in galaxies 
fainter than 0.06 L*^3 (---18 ABmag: see Figure E]). At 
~ — 18 AB mag, we find that galaxies have an average 
[/y-continuum slope /? of ~ — 2 to —2.5. 
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Fig. 10. — Bolometric luminosity versus observed UV luminosity Muy.ab- The top axis gives this UV luminosity in units of Lq. The 
blue points give the UV luminosities and equivalent bolometric luminosities (after dust correction) for individual sources in our z ^ 2.5 
(7-dropout and 2 ~ 4 B-dropout samples (see Figure [Sjl. The dust correction utilizes the observed (yV-continuum slopes /B and the Meurer 
et al. (l999) IRX-/3 relationship. Note that the red line is not a direct fit to the blue points (see §3.7 and FigureUJl. The red line shows the 
expected relation between bolometric luminosity and UV luminosity using the luminosity-dependent dust correction derived in §5.5 (see 
Eq. [2]l. These luminosity-dependent dust corrections are derived utilizing the correlation between (7l/-continuum slope 13 and Muv,AB 
we found for z ~ 2.5 (7-dropouts and z ^ A B-dropouts (see §3.9) and the Meurer et al. (1999) relationship. The hatched red region 
indicates the expected Itr scatter in the bolometric luminosities for a given UV luminosity (based upon the observed scatter of ~0.4 in the 
(7y-continuum slope at a fixed UV luminosity). Since it is unclear whether the Meurer et al. (1999) IRX-/3 relationship can be accurately 
used to estimate extinctions for galaxies with > lO'^^ Lq luminosities (e.g., Chapman et al. 2005; Reddy et al. 2006), caution should be 
used in interpreting the blue points or red lines which extend beyond lO'^^ Lq {green shaded region). The black dashed line shows the 
relationship between UV luminosity and /3 that one would derive utilizing the correlation Reddy et al. (2006) found between bolometric 
luminosity and dust obscuration (see §5.5). The shaded grey region corresponds to a ^ 0.4 dex scatter in the dust obscuration. Notice how 
the population of bluer, more dust-free dropout galaxies in our selections span a wide range in bolometric luminosity, extending all the way 
from 10^ Lq to 10^^ Lq where the LBG population transitions to a more dust-obscured ULIRG-like population. It is quite interesting that 
the UV luminosity where galaxies begin reaching bolometric luminosities of ~ 10^^ Lq - where dust obscuration becomes very important 
— is very close to L*_^ (Steidel et al. 1999: shown on the figure). The UV LF therefore begins to cut off at precisely the same luminosities 
as dust obscuration is becoming increasingly important in attenuating the light. This suggests that the value of L* at z ~ 2 — 4 may be 
set (in part) by the luminosities at which dust obscuration becomes dominant. It may also help to explain the mild evolution in the UV 
LF from 2 ~ 4 to z ~ 2 (see §5.5). 



These dust corrections are much more modest than 
have occasionally been assumed to calculate the total 
star formation rate density at early times, e.g., the fac- 
tor of 8 adopted by Giavalisco et al. (2004b) or the 
factor of ~3 adopted by Hopkins (2004) for their com- 
mon obscuration correction. Previously, Bouwens et al. 
(2006) and Lehnert et al. (2005) have argued that the 
dust extinction was likely less than a factor of < 2 at 
z ~ 5-6 and Reddy & Steidel (2009) argued that when 
integrated to lower luminosities the dust extinction was 
less than a factor of '--^1.5-2 at z^2-3 (adopting their 
Luminosity Dependent Reddening [LDR] model). The 
Reddy & Steidel (2009) LDR dust extinctions are some- 
what lower than what we estimate here, largely because 
they assume there is little dust extinction for galaxies 
with luminosities fainter than —18.5 (0.1 L*^^^). By 
contrast, the dust extinction we adopt come from the 
distribution of [/l^-continuum slopes we find and using 
Eq. [H While it is unclear which set of dust corrections 
is most accurate, it seems quite clear that the dust ex- 
tinction must be very low for galaxies with lower UV lu- 



minosities. This follows from two observational findings: 
(1) the much larger volume density (i.e., ~ 10~^ Mpc~^) 
of galaxies with low UV luminosities than ULIRG-type 
sources (i.e., where 10~'^'^ Mpc^^) and (2) the correlation 
between dust extinction and bolometric luminosity (e.g., 
Wang & Heckman 1996; Hopkins et al. 2001; Martin et 
al. 2005; Reddy et al. 2006; Buat et al. 2007; Burgarella 
et al. 2007). As a result of these two findings, we can 
conclude that essentially all of the galaxies with lower 
UV luminosities (>99%) do not have very high bolomet- 
ric luminosities, and therefore do not exhibit large dust 
extinctions. 

5.5. Bolometric Luminosities of Galaxies in our UV 
selections 

The viabihty of the Meurer et al. (1999) IRX-/3 rela- 
tionship (Eq. 1) we used to correct the observed UV 
luminosities for the dust extinction is a strong function 
of the bolometric luminosity of the sources. For exam- 
ple, Reddy et al. (2006) found that while the Meurer et 
al. (1999) relationship worked well between IO^^Lq and 
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10^^ Lq at z ~ 2, it failed badly for galaxies with lumi- 
nosities > lO^^L© (see also Chapman et al. 2005). 

To investigate the validity of the estimates we made 
in the previous section, it is interesting therefore to es- 
timate the bolometric luminosities for galaxies with spe- 
cific UV luminosities in our LBG selections. The results 
are shown in Figure [TOl We again take advantage of 
Eq. (1), as well as the relationships between the UV- 
continuum slope /? and UV luminosity (§3.9) to make 
this estimate. This results in the following relationship 
between the absolute magnitude of a galaxy in the UV 
and the bolometric luminosity Lf,oi at z ^ 2.5: 

Ltol = 10ll-67-0-58(A%v,A.,+21)^^ (2) 

The relationship at z ~ 4 is similar. We emphasize 
that this relationship does not generally hold for galaxies 
where Lboi > lO^^i© (e.g., Reddy et al. 2006; Chapman 
et al. 2005) or for galaxies which are particularly young 
(Reddy et al. 2006), because of difficulties in using the 
J/y-continuum slope /3 to correct for dust extinction in 
these cases. 

The above result is shown in Figure [10] with the red 
lines for the z ~ 2.5 and z ^ A samples. The hatched 
region shows the approximate scatter expected in the 
bolometric luminosity for a given UV luminosity as a re- 
sult of scatter in the t/V^-continuum slope f3. Also shown 
are the bolometric luminosities for individual dropout 
galaxies in the two samples, after correcting their UV 
luminosities for dust extinction using their measured /3's 
and Eq. (1). 

From FigurefTOland Eq.[2l it seems clear that we would 
expect a substantial fraction of our dropout sample to 
have bolometric luminosities less than 10^^ Lq. As such, 
we can plausibly make use of the IRX-/3 relationship to 
predict the mean dust extinction in galaxies - as we have 
done in the previous section. Of course, we might expect 
a few > lO'^^ Lq dust-obscured galaxies to be present 
in the normal dropout samples, with typical UV lumi- 
nosities and [/y-continuum slopes /? where Eq. [2] may 
not apply (e.g.. Chapman et al. 2005). However, given 
that the volume density of such galaxies is observed to 
be much lower (> 10-100 x) than that of typical dropout 
galaxies, we would expect Eq. [2] to be reasonably accu- 
rate in most cases. 

Examining Figure [TUl it is certainly quite striking that 
z ~ 2-4 LBGs start becoming extraordinarily rare bright- 
ward of ~ —21 or —22, at precisely the same luminosities 
(i.e., > 10^^ Lq) as we expect dust extinction to be- 
come increasingly important in galaxies. This suggests 
the abrupt cut-off in the UV LF at z^2-4, i.e., L^y, is 
set by dust extinction - or more precisely the luminos- 
ity at which dust obscuration becomes so significant as 
to offset the increase in energy output from stars. We 
can examine this by using the observed correlation be- 
tween dust extinction and bolometric luminosity in the 
blue star-forming population at z ~ 2 (Reddy et al. 2006) 
log Lboi = (0.62 ± 0.06) log -I- (10.95 ± 0.07) and ap- 
plying it to Figure [TOl Adopting the observed scatter of 
A log ~ 0.4, we show the envelope of Muv and Lf,oi 
values expected for a wide range in Lboi- It seems clear 
from this exercise that we would not expect galaxies to 
have UV luminosities much brighter than —22 (~2 L^y) 
and that dust effectively sets the cut-off in the UV LF. 



This may explain why the value of L^y (e.g., Steidel et 
al. 1999; Reddy et al. 2008) does not evolve much from 
z ^ A to z ^ 2. 

5.6. High-Redshift Galaxies from the Literature that 
would seem to show Substantial Dust Extinction 

The dust extinctions we infer are substantially less 
than have been inferred in some studies of specific high 
redshift galaxies. One such example is the z — 6.56 
galaxy behind Abell 370 (HCM6A: Hu et al. 2002). This 
source was found to have a significantly larger flux in the 
4.5^ band than in the 3.6/i band - which may be due 
to substantial Ha emission falling within the 4.5/i band 
(Chary et al. 2005). If Ha is in fact the explanation, the 
inferred star formation rate for this source from this Ha 
emission is ^ 10 x higher than would be inferred from the 
UV or Lya emission alone. Detailed stellar population 
modelling on the optical, near-IR, and IRAC imaging 
data of IICM6A suggest an extinction factor of ^ 11 x 
{Auv — 2.6 mag) which is much larger than what we 
infer on average for z ~ 6 i-dropouts (Chary et al. 2005; 
Schaerer & Pello 2005; see Table [5]). However, in subse- 
quent and much deeper IRAC observations of HCM6A 
(K. Lai et al. 2009, in prep), no significant difference be- 
tween the fluxes of HCM6A in the 3.6/i and 4.5/i bands 
is found, leaving us with no reason to believe this source 
shows significant Ha emission and hence no reason to 
suspect that its UV fight is significantly obscured. Direct 
observations of this source with the MAMBO-2 bolomet- 
ric array at 1.2mm also set useful upper limits on its dust 
obscuration (Boone et al. 2007). 

Other examples might include the many high redshift 
QSOs discovered by the Sloan Digital Sky Survey (SDSS: 
Fan et al. 2006). The inferred dust masses of SDSS QSOs 
can be quite substantial (e.g., J1148-I-5251: Bertoldi et 
al. 2003), providing clear evidence that specific high red- 
shift sources can show modest to sizeable dust extinction. 
Nonetheless, we would not expect QSOs like J1148-I-5251 
to be representative of z ^ 6 galaxies, associated as they 
are with very massive halos and thus in very rich en- 
vironments. Galaxies in such environments are known 
to show much greater evolution than galaxies in more 
typical environments (e.g., Tanaka et al. 2005) and show 
much greater metal enrichment. 

6. POSSIBLE IMPLICATIONS FOR THE STAR 
FORMATION HISTORY 

In §5, we used the observed distribution of UV- 
continuum slopes (3 to estimate the dust extinction in 
star- forming galaxies at z ~ 2 — 6. In the current sec- 
tion, we apply these dust extinctions to UV LFs derived 
at z '--^ 2 — 6 to estimate the SFR density from typical 
star-forming galaxies (§6.1). In §6.2, we use recent search 
results for ULIRGs to estimate their likely contribuion to 
SFR density (since they will not be accurately accounted 
for by the dust-corrected LBG results). In §6.3, we use 
these two results to estimate the bolometric flux from 
SFR to be output in the UV and IR. Finally, in §6.4, we 
put together these results to estimate the SFR density 
from all star-forming galaxies at z ~ 2 — 6. 

6.1. The Inferred Star Formation Rate Density at 
z^2-6 
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TABLE 6 

Inferred Star Formation Rate Densities from 





Z 


'2-6 LBG SAMPLES. ''''^ 


Dropout 




logiQ SFR density (Mq Mpc-3 yr^^) 


Sample 


< z > 


L > 0.3L*^3 


L > 0.04L*^3 








Uncorrected 


u 


2.5 


-1.70±0.03 


-1.36±0.03 


B 


3.8 


-1.81±0.05 


-1.48±0.05 


V 


5.0 


-2.15±0.06 


-1.78±0.06 


i 


5.9 


-2.31±0.08 


-1.83±0.08 


z 


7.4 


-2.58±0.21 








Dust-Corrected 


U 


2.5^^ 


-0.93±0.03 


-0.78±0.03 


B 


3.8 


-1.05±0.05 


-0.90±0.05 


V 


5.0 


-1.72±0.06 


-1.57±0.06 


i 


5.9 


-2.10±0.08 


-1.70±0.08 


z 


7.4 


-2.37±0.21 





^ Based upon LF parameters in Table 2 of Reddy & Steidel 
(2009), Table 7 of Bouwens et al. (2007:see §3.5), and Table 
4 of Bouwens et al. (2008). 

^ The SFR density here tabulated in terms of the Salpeter 
IMF. Expressing these results in terms of a Kroupa (2001) 
IMF, one should divide the results given here by a factor of 
~1.7. 

" We adopt the Reddy & Steidel (2009) UV LF at z = 3 for 
computing the SFR density at z ^ 2.5. 

It makes sense for us to use the effective dust extinction 
inferred in the previous section (Table O to revise our 
estimates the star formation rate density at z ~ 2-6. As 
per our calculation above, we adopt the UV LFs derived 
by Reddy & Steidel (2009) at z ~ 3 and Bouwens et 
al. (2007) in calculating the SFR densities. In order to 
convert luminosities in the UV into estimates of the SFR 
we use the transformation presented by Madau et al. 
(1998: see also Kennicutt 1998): 



Luv — 



/ SFR 



8.0 X lO^Vgss^^Hz"^ 



(3) 



where a 0.1-125 M0 Salpeter IMF and a constant star 
formation rate of > 100 Myr are assumed. Conversion 
to a Kroupa (2001) IMF would result in factor of ^1.7 
smaller SFR estimates. 

In view of the young ages (~ 10-50 Myr) inferred for 
many star-forming galaxies at z~5-6 (e.g., Yan et al. 
2006; Eyles et al. 2007; Verma et al. 2007), it is clear 
that adopting a simple formula like Eq. [3] (which require 
the assumption of > 100 Myr ages) would result in a 
systematic underestimate of the SFR density of the uni- 
verse at very early times (e.g., Verma et al. 2007; Reddy 
& Steidel 2009; Stanway et al. 2005). In any case, given 
our uncertain knowledge of how the age distribution of 
star-forming galaxies varies as a function of redshift and 
luminosity, it is difficult to know precisely how much we 
must revise our estimates upward. Working with a sam- 
ple of bright z ~ 5 galaxies, Verma et al. (2007) estimate 
that the actual SFRs may be twice as large as what Eq. [3] 
would suggest. Of course, this effect may be partially 
offset by the fact that for these same young star-forming 
systems that standard dust corrections (e.g., following 
Meurer et al. 1999) seem to overpredict the SFRs by fac- 
tors of >2 (Reddy et al. 2006; Siana et al. 2008). In any 
case, going forward, we will need to obtain tighter con- 
straints on the age distribution of galaxies as a function 
of redshift and especially to very low luminosities (since 
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Fig. 11. — The star formation history inferred using the dust 
corrections determined here (upper set of points: orange shaded 
region) and without those dust corrections (lower set of points: 
blue shaded region). These dust corrections (Table O are applied 
to the luminosity density determinations at 2 > 2 by Reddy & 
Steidel (2009), Bouwens et al. (2007), and Bouwens et al. (2008). 
Otherwise, the dust corrections of Schiminovich et al. (2005) are 
assumed. The top panel is based upon the total UV luminos- 
ity density integrated to 0.3 L*_^ and the bottom panel is based 
upon the total UV luminosity density integrated to 0.04L*_.^. The 
rest-frame UV continuum luminosity density is converted to a star 
formation rate density assuming a constant > 10* yr star formation 
model and a Salpeter (1955) IMF (Madau et al. 1998). Conver- 
sion to a Kroupa (2001) IMF would result in a factor of ~1.7 (0.23 
dox) decrease in the SFR density estimates given here. Assum- 
ing younger ages for star-forming galaxies in high-redshift sam- 
ples would increase the SFR density estimates by a similar factor 
(Verma et al. 2007). Shown are the luminosity density determi- 
nations by Schiminovich et al. (2005: black hexagons), Reddy & 
Steidel (2009: green crosses), Bouwens et al. (2007: blue/red cir- 
cles), and Bouwens et al. (2008: blue/red circles and la upper 
limit). The open symbols show the star formation rate densities 
that would be inferred if we applied the same dust correction to 
> 0.04L*_3 as we apply to > 0.3Z;*_3. It is clear then that ignor- 
ing the fact that lower luminosity galaxies have much lower values 
for the dust extinction will result in an overestimate of the SFR 
density. See also Table [5] Table 5 and Figure 10 from Reddy & 
Steidel al. (2009) make a similar point. The equivalent star forma- 
tion rate density (UV luminosity density) is also shown assuming 
no extinction correction (lower set of points). 

faint galaxies provide the dominant contribution to the 
SFR density). 

Our results are presented in Figure [Tl] and Table [6l 
These dust-corrected SFR densities are somewhat lower 
than presented by Bouwens et al. (2008). This is be- 
cause the average dust extinction including the lowest 
luminosity galaxies is less than when only considering 
L* galaxies. We will discuss these results in more depth 
in §6.4 after including the contribution from ULIRGs. 

6.2. Contribution of ULIRGs to the Star Formation 
Rate Density at z ^ 2-6 

Of course, we would expect the extinction estimates 
above to underestimate the effective extinction, given 
that Eq.[T]does not accurately describe galaxies with very 
high SFRs, dust extinctions, and bolometric luminosi- 
ties (e.g., Elbaz et al. 2007), particularly those systems 
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Fig. 12. — Fraction of the total luminosity density from stars 
at 2 ~ 2-6 emitted by galaxies in the rest-frame UV and in the 
mid-IR/far-IR (after dust reprocessing). In calculating the energy 
output in the UV vs. IR, we adopt the extinction estimates in Ta- 
bleland assume that the total energy output from ULIRGs occurs 
in the infrared (see §6.3). A dotted line is included to show the 
relative contribution of LBGs and ULIRGs to the mid-IR/far-IR 
energy output. The contribution of star-forming galaxies is in- 
cluded all the way down to zero luminosity. Most of the energy 
output at 2 ~ 2-4 would appear to occur in the IR, while at 2 > 5 
the energy output occurs mostly in the UV. Note that despite 
most of the energy output at z ~ 2.5 — 4 occurring at IR wave- 
lengths, the SFR densities we derive from LBG surveys (after dust 
corrections ) sh ould account for > 80% of the SFR density at z > 3 
(see Figure [Ml . 

where Li,oi > 10^^ Lq. The SFR from such galaxies can 
be better accounted for by a direct consideration of IR 
LFs (e.g., Caputi et al. 2007).'* Reddy & Steidcl (2009) 
have aheady performed such a calculation at z ^ 2, in- 
tegrating the z ~ 2 IR LF of Caputi et al. (2007) down 
to bolometric luminosities of lO-'^^L©.^ Reddy & Stei- 
del (2009) infer a SFR density of 0.03 Afoyr-iMpc"^ 
(Kroupa IMF) at z ~ 2 from this population. Obviously, 
there are some modest uncertainties as to the exact bolo- 
metric luminosity down to which one should integrate 
the Caputi et al. (2007) LF (to correct for the ULIRG 
population not reasonably accounted for by the Meurer 

^ While there has been some suggestion that galaxies whose IR 
luminosities exceed that predicted by the IRX-;9 relationship fail 
due to a contribution from an obscured AGN (e.g., Daddi et al. 
2007) and therefore the SFR derived from the IRX-/3 is correct 
after all, it has been found (e.g.. Murphy et al. 2009) that this is 
only partially the explanation, and that part of the IR excess would 
appear to be due to star formation. We therefore consider it safer 
to include the contribution from > 10^^ L© IR-luminous galaxies 
explicitly through the use of IR surveys. 

^ There has been sizeable variation in the literature in the rela- 
tionships used to convert from 24/j fluxes to the bolometric flux in 
the IR, (e.g., Bavouzet et al. 2008 vs. that implied by the Chary 
& Elbaz 2001 templates found to work in the local universe), and 
therefore large differences in the reported LFs in the IR (e.g., Perez- 
Gonzalez et al. 2005 vs. Caputi et al. 2007). Fortunately, the 
situation seems to have become more clear from the observations, 
and Caputi et al. (2007) adopt an apparently more justified linear 
relationship (Bavouzet et al. 2008) to convert between the 24/im 
flux and the bolometric IR flux. 



et al. 1999 relation), and so this SFR density could be 
larger. For example, assuming we integrated the Caputi 
et al. 2007 LF down to 10"-® Ltoi instead, the ULIRG 
contribution could increase by ~50%. 

Despite the existence of these possible uncertainties, 
we adopt as our baseline values at z ~ 2.5 the contri- 
bution Reddy & Stcidel (2009) derived from the z 2 
Caputi et a. (2007) IR LF. We therefore correct our SFR 
density estimates at 2 ^ 2.5 in Table [6] for this popula- 
tion. The results are given in the lower section of Table[7l 
The obscured population only has a modest effect on the 
dust extinction (and SFR density), increasing it by ~20% 
(see also discussion in Reddy & Steidel 2009). 

Accounting for the effect of such ULIRGs on the SFR 
density (and extinction estimates) at z > 4 is a much 
less certain enterprise. Not only is it difficult to probe 
to even 10^^ Lq at z A with current instrumentation, 
but conducting follow-up observations on specific z > A 
sources (e.g., to determine redshifts) is also quite chal- 
lenging at present. As a result, only a handful of high 
quality z > 4 ULIRG candidates (e.g., Capak et al. 2008; 
Daddi et al. 2009; Wang et al. 2009; Schinnerer et al. 
2008; Dunlop et al. 2004) currently exist, and essentially 
all of them are at z ^ 4. Perhaps the deepest, most well- 
defined selection of such sources is found in the deep 
^50-100 arcmin^ submm SCUBA field lying within the 
HDF-North GOODS area (Pope et al. 2006). We will 
use those z ~ 4 sources to set a lower limit on the con- 
tribution of highly obscured galaxies to the SFR density 
at z ^ A. Starting with the ^5-8 ULIRGs within this 
field that are plausibly at z ^ 4 (e.g., Daddi et al. 2009; 
Mancini et al. 2009) and then dividing by the comoving 
volume at z ~ 4, one derives a SFR density of ^0.01-0.02 
M0yr~*Mp~'^ (depending upon what fraction of the can- 
didates are actually at z ~ 4). Again, we can correct our 
z ^ 4 SFR density estimates for this population. Not 
surprisingly, this population only contributes 10% to 
the SFR density (see Figure [Ti)). 

Of course, we might expect there to be highly obscured 
galaxies below the 10*^ Lq luminosity limits that also 
contribute substantially to the SFR density relevant to 
these probes, so one might argue that 10% is just a lower 
limit. However, as we have already argued, we should 
be able to account for the star formation of these lower 
luminosity systems from LBG selections. The results of 
Reddy et al. (2006) suggest that the f/lZ-corrected SFRs 
of galaxies can be accurately estimated at z ~ 2 to 10^^ 
Lq, and we would expect these corrections to be even 
easier to perform at z '--^ 4 than at z ~ 2 - given the 
trend of high redshift galaxies to show less and less dust 
extinction at a fixed SFR at increasingly high redshifts 
(e.g., Reddy et al. 2006; Buat et al. 2007; Burgarella et 
al. 2007). 

Given that > 10*^ Lq ULIRGs only contribute --20% 
of the SFR density at z '--^ 2-3 and perhaps ^^10% at z ~ 4, 
one might expect the contribution of such a population 
at z ~ 5-6 also to be quite small. We can obtain a crude 
estimate of this fraction, if we assume that IR luminous 
galaxies have a bolometric luminosity function similar to 
that derived at z ~2 (e.g., by Caputi et al. 2007 where 
-^boi ~ 10*^ ®^°'*) and then evolve in a very similar way 
to the UV LF at z > 3. We therefore assume that 
L'^^j scales with redshift in the same way as L^^jy, i.e.. 
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TABLE 7 

Inferred Star Formation Rate Densities, including 

THE contributions FROM HIGHLY DUST OBSCURED GALAXIES 

(§6.2-6.3).^-'^ 



10 54 3 2 



t (Gyr) 

1 0.6 



0.4 



Dropout 



loglO SFR density (Mq Mpc"^ yr"^) 



Sample 


< z> 


L > 0.3L*^3 


L > 0.04L*^3 






Dust-Corrected 


U 


2.5'= 


-0.93±0.03 


-0.78±0.03 


B 


3.8 


-1.05±0.05 


-0.90±0.05 


V 


5.0 


-1.72±0.06 


-1.57±0.06 


i 


5.9 


-2.10±0.08 


-1.70±0.08 


z 


7.4 


-2.37±0.21 








Dust-Corrected -|- ULIRC^ 


U 


2.5'= 


-0.77±0.03 


-0.66±0.03 


B 


3.8 


-0.98±0.05 


-0.85±0.05 


V 


5.0 


-1.66±0.06 


-1.53±0.06 


i 


5.9 


-2.07±0.08 


-1.69±0.08 


z 


7.4 


-2.37±0.21 





Based upon LF parameters in Table 2 of Reddy & Steidel 
(2009), Table 7 of Bouwens et al. (2007:see §5.6), and Table 
4 of Bouwens et al. (2008). 

The SFR density here tabulated in terms of the Salpeter 
IMF. Expressing these results in terms of a Kroupa (2001) 
IMF, one should divide the results given here by a factor of 
~1.7. 

" We adopt the Reddy & Steidel (2009) UV LF at z = 3 for 
computing the SFR density at z ^ 2.5. 

At 2 ~ 2.5, we include the SFR density contribution from 
> 10^^ Lq ultraluminous IR galaxies (ULIRG: §6.2) by using 
the 2 ~ 2 IR LF of Caputi et al. (2007) as Reddy & Steidel 
(2009) do. At 2 ~ 4, we estimate their contribution based 
upon a small sample of 2 ~ 4 ultra-luminous IR galaxy candi- 
dates (see Daddi et al. 2009) within the HDF-North GOODS 
SCUBA supermap (Pope et al. 2006). Based upon the 2 ~ 4 
result, we assume that ULIRGs provide only a ~5% contribu- 
tion to the SFR density at 2 ~ 5-6. We believe this assump- 
tion is a reasonable one given the trends towards much bluer 
l/V-continuum slopes at high redshifts seen in both LBG and 
BBG samples (§4.2) and thus apparently smaller fraction of 
galaxies with substantial dust obscuration. 

as io-0-4(0-36(2-3)) _ (1/1.39)^-3 to match the observed 
M^y{z) = Ml}y{z = i) + OM{z - 4) evolution in the 
UV LF from z - 4 to z - 7 (e.g., Bouwens et al. 2007; 
Bouwens et al. 2008).® We then compare the fractional 
SFR density in LhoZ > IO^^Lq galaxies with that from 
the star-forming population as a whole. We find that the 
fractional SFR density in ULIRGs is 6% at z ~ 5 and 
2% at z 6. For these assumptions, the fractional SFR 
density in Lhoi > IO^^Lq galaxies is 23% and 10% at 
z ~ 3 and z ^ 4, respectively, which would seem to be 
a good match to the observations. We emphasize that 
this very low ULIRG estimate at z >5 is consistent with 
the lack of very red J/V'-continuum slopes /3 at z ~ 5-6 
(Figure [S] which is in contrast to the tail towards redder 

[i.e., 1 to 0] /3's at z - 4). 

Recent work on GRBs provide strong support for the 
general conclusion that highly obscured ultra-luminous 
IR galaxies contribute at most a small fraction of the 

^ Note that here we are implicitly assuming that L^^j evolves 
from 2 ~ 4 to 2 ~ 3 despite the lack of evolution in L'^y over 
this redshift range (e.g., Reddy et al. 2008). Such a brightening in 
L* would be a natural continuation of the evolution in the UV LF 
from 2 ~ 7 to 2 ~ 4 (e.g., Bouwens et al. 2008). The reason we 
may not observe this brightening in the UV LF (i.e., L^y) beyond 
2 ~ 4 is because of the increasing importance of dust extinction 
in the most bolometrically luminous galaxies (e.g., see §5.5 where 
we argue that dust extinction may be partially responsible for the 
cut-off in the C/V LF at 2 ~ 2 - 4). 
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Fig. 13. — Star formation history inferred with (upper set of 
points, darker orange contours) and without (lower set of points, 
lighter orange contours) a contribution from highly obscured ultra- 
luminous infrared bright (> IQ^^Lq) galaxies. Otherwise similar 
to Figure nn The top panel is based upon the total UV luminos- 
ity density integrated to 0.3 L*_„ and the bottom panel is based 



upon the total UV luminosity density integrated to 0.04L* 



At 



2 ~ 2.5, we follow Reddy & Steidel (2009) in using the 2 ~ 2 IR LF 
of Caputi et al. (2007) to include the contribution from these ultra 
luminous IR galaxies, while at 2 ~ 4, we estimate their contribu- 
tion based upon ~5-8 good 2 --^ 4 ULIRG candidates (see Daddi et 
al. 2009) within the HDF-North GOODS SCUBA supermap (Pope 
et al. 2006). See also Table [7] and §6.2. The darker orange con- 
tours at 2 ~ — 1.2 show the SFR density derived using the deepest 
mid-IR/far-IR observations over the GOODS and Far Infared Ex- 
tragalactic Legacy (FIDEL) fields (Magnelli et al. 2009). While 
IR bright galaxies appear to add significantly to the SFR, density 
at low redshift in the upper panel, this is because the lower set of 
points/orange contours only include galaxies at the bright end of 
the UV LF (corresponding to the O.SL*^,^ limit) and the UV LF 
there cuts off at much fainter magnitudes than at higher redshifts 
(compare this situation to the lower panel where the contribution 
of the mid-IR/far-IR galaxies is much less on a percentage basis). 
From this figure, it is clear that the contribution of the luminous 
IR sources to the total s tar formation rate density at 2 > 2 is only 
modest (see also Figure [T4)l . 

SFR density at z > 4. Since GRBs are perhaps a good 
tracer of massive star formation in the high redshift uni- 
verse (e.g., Yiiksel et al. 2008; Li 2008; Kistler et al. 
2009) and may not be biased to lower metallicity envi- 
ronments (Prochaska et al. 2007), we would expect GRBs 
to occur in galaxies roughly in proportion to their share 
of the star formation activity in the high redshift uni- 
verse. Follow-up work done on GRB host galaxies show 
that these systems are quite faint at both rest-frame UV 
(Chen et al. 2009) and optical wavelengths (Chary et 
al. 2007a). We can therefore conclude that most of the 
star formation at high redshift occurs in galaxies that are 
similarly faint. 

6.3. Fractional Light Output in the UV and 
mid-IR /far-IR 
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Fig. 14. — Fraction of the total SFR density at high redshift 
in lower luminosity LBGs (< 0.3-L*_3), higher luminosity LBGs 
(> 0.3L*_3), and ultra-luminous (> IO^^Lq) infrared galaxies 
(ULIRGs). While there will obviously be some overlap between 
LBG and ULIRG selections, we consider the ULIRG population 
here separately since mid-IR/far-IR data are required to fully ac- 
count for the SFR in these galaxies. The contribution of star- 
forming LBGs is included all the way down to zero luminosity. 
These fractions are derived based upon the discussion i n §5 . 4, § 6.1, 
and §6.2, and estimates in Tables K\ and [7l and Figures 111! and [131 
It is clear that virtually all the SFR density at z ^3-4 can be 
accounted for from LBG selections. 

Before using the extinction estimates for LBG samples 
(Tabic m §5.4) and the aforementioned estimates of the 
tlLIRG contribution (§6.2) to estimate the SFR density, 
we first use them to assess the relative flux output from 
galaxies in the UV and that output in the mid-IR/far- 
IR, after dust reprocessing. This entire topic is impor- 
tant in thinking about how best to nail down the star 
formation rate density and energy output in the high- 
redshift universe and for our thinking about the contri- 
bution of the high-rcdshift universe to the extragalactic 
background light. 

In calculating the relative flux output in the UV and 
mid-IR/far-IR, we assume that essentially all the light 
from the ULIRG population is emitted in the mid-IR/far- 
IR and adopting the extinction estimates from Table [5] 
for the LBG population. We include a plot of this frac- 
tion versus redshift in Figure 1121 Not surprisingly the 
majority of the light output from star-forming galaxies at 
z ~ 2.5 — 4 occurs in the mid-IR/far-IR (after dust repro- 
cessing) as has been emphasized repeatedly (e.g., Hughes 
et al. 1998; Perez-Gonzalez et al. 2005). However, as we 
will note in §6.4, it appears that we can account for most 
of this energy output using the LBG population. Again, 
the situation appears to be quite different at z > 5, with 
most of the light output occurring in the UV. 

6.4. The Star Formation Rate Density at z ^ 2-6 
Including ULIRGs 

Including the contribution from highly obscured ultra- 
luminous IR galaxies (§6.2) in with that derived from the 
UV LFs (Tableland Figure fTTj) . our best estimates for 



the star formation rate densities are given in Table [7] and 
Figure [131 The ULIRG population increases the SFR 
density estimated at z ~ 2.5 and z ~ 4 (to faint limits) 
by ~ 20% and 10%, respectively. We show the fractional 
contribution of ULIRGs to the SFR density at z ~ 2-6 in 
Figure [H 

The present dust-corrected SFR densities are less than 
has been considered in many previous studies. This is 
largely because many previous studies corrected the star 
formation rate density by factors of ~ 5 — 8x (e.g., Gi- 
avalisco et al. 2004b) - which is the dust correction rele- 
vant for L* galaxies. This has important implications for 
studies (e.g., Hopkins & Beacom 2006) featuring compar- 
isons of the stellar mass density with the star formation 
rate density. 

The claim has been that the star formation rate den- 
sity - after integration - exceeds the stellar mass density 
by factors of '--^ 3 — 4 (e.g., by Hopkins & Beacom 2006; 
Wilkins et al. 2008). Adopting a small dust correction (as 
is appropriate including the effect of the lower luminos- 
ity galaxies) substantially reduces the inferred star for- 
mation rate density at high redshift (by factors of ~ 1.1, 
~1.1, ~2.1, and '--^2.5 at z ~ 2.5, z ~ 4, z ~ 5, and 
z ~ 6, respectively relative to that used by Wilkins et 
al. 2008). This should help somewhat to resolve the ap- 
parent discrepancy in comparisons with the stellar mass 
density inferred at high redshift (see also discussion in 
Reddy & Steidel 2009). Another consideration that may 
help in this regard is the stellar mass in lower luminosity 
galaxies. Reddy & Steidel (2009) argue that a considera- 
tion of the SFR density and stellar mass densities to the 
same limiting luminosities largely resolves the apparent 
discrepancy between the inferred stellar mass density and 
that obtained by integrating SFR history (approximately 
doubling the inferred stellar mass densities at z ~ 2-3). 

What stands out in the above discussion is the role of 
lower luminosity galaxies and their apparently sizeable 
contribution to various volume-averaged quantities like 
the SFR density and stellar mass density. Clearly, any 
reasonable examination of these quantities (or compar- 
isons between the stellar mass density and the integral of 
the SFR density) must accurately account for this popu- 
lation or risk ignoring those galaxies that are most central 
to making sense of the cosmic evolution. 

7. SUMMARY 

We use the deep optical and near-IR imaging data over 
the HUDF and other deep, wide-area fields to quantify 
the distribution of t/F-continuum slopes (3 (i.e., f\ (x \^) 
for star-forming galaxies over a wide range in redshift 
(z'-2-6) and luminosity (O.lLJ^g to 2Ll^^). z-2-6 
galaxies are selected through a U, B, V, and i dropout 
technique while J/y-continuum slopes (1600-2300A rest- 
frame) are derived from the optical and near-IR broad- 
band colors. We then corrected the distribution of UV- 
continuum slopes for observational selection and photo- 
metric errors and tabulate these slopes as a function of 
both redshift and luminosity (Tableland Figure [3]). We 
then discuss possible interpretations of the trends we find 
(§4). 

Our conclusions are as follows: 

• The C/y-continuum slope distribution of UV bright 
L*^3 galaxies over the range z ~ 2-6 has a mean /] 



ranging from —1.2 to —2.4, with a dispersion of 
^--^0.4 (Tableland Figure [S]). As found in previous 
studies, we find that the mean C/T^-continuum slope 
(3 is bluer (by ^ 0.5) at z ^ 6 than it is at z ~ 3 — 4. 
We also find that the mean UV-continuum slope 
(3 is bluer (by ^0.5) at lower luminosities than it 
is at higher luminosities (see also Meurer et al. 
1999). In doing so, we establish the following 
correlation between (3 and Muv ab (§3-9): (3 = 
(-0.20±0.04)(M[/yAB + 21)-(l'.40±0.07±0.15) 
at z - 2.5 and/3 = ('-0.15 ± 0.01)(Mc/y,AB + 21) - 
(1.48 ± 0.02 ± 0.15) at z ~ 4. The present quantifi- 
cation of these distributions is essential for accurate 
LF determinations at high redshift (see §4.3). 

The dropout color selections identify star-forming 
galaxies at z ~ 2-6 with t/F-continuum slopes (3 as 
red as (3 ~ 0.5 (see Figures [2] and |4]) . By contrast, 
the distribution of f/y-continuum slopes (3 that we 
observe has a mean value of ~ —1.2 to —2.4 with 
a 1(7 dispersion of 0.4 - bluer than the color se- 
lection limit (i.e., /3 ~ 0.0 - 0.5) by A/3 > 1.5. 
If the distribution of [/^-continuum slopes has a 
simple form (i.e., has only a single mode), this sug- 
gests that our high-redshift dropout selections are 
largely complete (see also §4.1 in Bouwens et al. 
2007), and there is not a large population of dust- 
reddened galaxies at z > 4. This deficit of red 
galaxies is particularly prominent relative to the 
selection function at z > 5 (see Figure [6] §4.2). 
Independent evidence for this come from Balmer 
break selections at z ~ 3 — 4.5 (which are biased 
towards older and more massive galaxies) where 
the distribution of t/y-continuum slopes /3 is also 
found to be quite blue (Brammer & van Dokkum 
2007: see §4.2). 

Stellar population models were used to investigate 
the effects that changes in dust, age, metallicity, 
and IMF would have on the [/^-continuum slope 
(3. We found that a factor of ~2 (0.3 dex) decreases 
in the dust, age, and metallicity would make the 
f/y-continuum slope (3 bluer by ^0.35, ~0.1, and 
~0.05, respectively. Changing the slope of the IMF 
of galaxies by ~ 0.5 from Salpeter has a minimal 
effect, changing /3 by ^0.1. The UV -coTiiumnm 
slope /3 appears to be most sensitive to changes in 
the dust extinction. Dust is likely the most signifi- 
cant cause for changes in (3 as a function of redshift 
and luminosity (%^.^). 

Assuming that the observed correlation between 
dust extinction and C/y-continuum slope /3 at z ~ 
0-2 (e.g., Meurer et al. 1995; Meurer et al. 1999; 
Burgarella et al. 2005; Laird et al. 2005; Reddy et 
al. 2006; Dale et al. 2007) holds out to z ~ 6, we es- 
timate the probable dust extinction in z ~ 2-6 LBG 
selections (see Table [5]). For > 0.3i*^3 galaxies at 
z ~ 2.5, our estimate is 6.0±2.5 - very similar to 
the values derived by Reddy et al. (2006) and Erb 
et al. (2006b) using a variety of different multiwave- 
length data and SFR calibrators, but ~3x higher 
than CariUi et al. (2008) estimated for >0.2L*^3 
z ^ 3 {/-dropouts by stacking radio observations 
over the COSMOS field. At zr^5-6 and to much 



lower luminosities (> 0.04^*^3^, the inferred dust 
extinction is much less, i.e., <2-3x. 

• Since lower luminosity galaxies dominate the lumi- 
nosity density in the UV at z > 2 (Bouwens et al. 
2006, 2007; Yan & Windhorst 2004; Beckwith et 
al. 2006; Reddy & Steidel 2009: see Figure E]), es- 
tablishing the dust extinction of lower luminosity 
galaxies is crucial for estimates of the total star 
formation rate density. As a result of the very 
low dust extinctions inferred for lower luminosity 
galaxies, we expect the average dust extinction for 
the star-forming population at z ^ 2-6 integrated to 
very low luminosities to be quite small (< 2x: see 
Table [5^. This is similar to what Reddy & Steidel 
(2009) concluded from the observations and vari- 
ous physical arguments (see Table 5 and Figure 10 
from that work). 

• We establish a conversion between UV lumi- 
nosities and bolometric luminosity based on the 
aforementioned correlation between dust extinc- 
tion and (3 and the observed correlation between 
(3 and Muv,ab: Lboi = lO^^-e7-o.58{Muv,AB+2i) 
(though we emphasize that this relationship may 
not work for galaxies where Li,oi > 10^^ -^o)- We 
remark that it is striking that the UV LF cuts off at 
luminosities (~l-2 L^y) that approximately cor- 
respond to that of ULIRGs (^ 10^^ Lq: see Fig- 
ure [10]). This suggests that L^jy at z Z is set by 
dust extinction - or more precisely the luminosity 
at which dust obscuration becomes so significant as 
to offset the increase in energy output from stars. 
This may explain why the value of L^y (e.g., Stei- 
del et al. 1999; Reddy et al. 2008) does not evolve 
much from z ^ 4 to z ^ 2 (see §5.5). 

• We correct our SFR density estimates upperward 
to account for the contribution from dust obscured 
ULIRGs. At z ^ 2.5, we follow the analysis of 
Reddy & Steidel 2009 in using the z - 2 IR LF of 
Caputi et al. (2007) to account for the contribution 
from these galaxies and make a similar correction 
at z ~ 4 based upon ~ 5 — 8 plausible submm 
candidates at z ~ 4 in the HDF-North GOODS 
SCUBA field discussed by Daddi et al. (2009: see 
§6.2). We find that this obscured population only 
increases the SFR density (and dust extinction) es- 
timates by ^20% at z ^ 2.5 and by just ^ 10% at 
z ^ 4 (Figure \14^ . Given the evolution towards 
bluer UV-continuum slopes (3 at high redshifts in 
both LBG and BBG selections (^4-^)j W6 argue that 
the contribution of dust-obscured ULIRGs to the 
SFR density at z > 5 is even less than at z ^ 4 
(i.e., < 10%). The model we construct in ^6.2 (as- 
suming that L^^i evolves similarly to L^y) suggests 
that the fractional contribution of ULIRGs to the 
SFR density is only 6% and 2% at z ^ b and z ~ 6, 
respectively. 

The distribution of C/y-continuum slopes at z ~ 2-6 pro- 
vides us with a key window of information for study- 
ing the formation and evolution of star-forming galaxies 
from very early times. The very blue colors of galaxies 
at z ~ 5-6 and at lower luminosities suggest that they 
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are indeed much less evolved than their lower redshift 
or higher luminosity counterparts. These results suggest 
that the dust corrections at high redshifts are not par- 
ticularly large and that LBG (dropout) selections should 
give a substantially complete census of the SFR density 
at z > 5 and probably even at z > 4. Moreover, highly 
dusty galaxies seem unlikely to provide a significant con- 
tribution to the SFR density at these early times. 

In the future, we expect to be able to constrain the 
distribution of C/V-continuum slopes p at z ~ 1 — 2 
and z ^ 5 — 8 with much more precision and accuracy 
than is possible at present. These advances should come 
using data from the HST WFC3/UV and WFC3/IR 
instrument to be installed on HST. This instrument will 
allow us to collect deep UV data needed to efficiently 
select star-forming galaxies at z ~ 1 — 3 and the near-IR 
data needed to estimate the t/T^-continuum slopes (3 for 
z ~ 4 — 8 galaxies. 
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Fig. 15. — Expected W colors for star-forming galaxies in our z ~ 2.5 [/, 2 ~ 4 B, z ~ 5 V, and z ~ 6 i dropout selections. The black 
lines (from bluest to reddest) correspond to power-law SEDs with /3's of —2.2, —1.5, and —0.8, respectively. The red lines (from bluest to 
reddest) correspond to SEDs computed assuming an e"*/"^ star formation history, t = 10 Myr, r = 70 Myr, [Z/Zq] = —0.7, the Bruzual 
& Chariot (2003) libraries, a Salpeter IMF, and dust extinctions E{B - V) of 0.0, 0.15, and 0.3, respectively (Calzetti et al. 2000). These 
latter models have effective C/V-continuum slopes /3 of —2.2, —1.5, and —0.8 over the range I6OOA to 23OOA - the same as the black 
lines. However, the red curves should be much more realistic - being based upon stellar population models (Bruzual &; Chariot 2003). We 
reference the measured colors from each of our dropout samples to this fiducial set of SEDs (red lines) to account for intrinsic differences 
in the rest-frame wavelengths of the broadband imaging data available for each of our samples (Table [2] and §3.4). Also shown (dashed 
and dotted lines, respectively) are the predicted colors for SEDs computed assuming a young stellar population (a 30 Myr constant star 
formation rate model) and a somewhat older stellar population (a model for which star formation proceeded at a constant rate for 100 
Myr and then ceased for 10 Myr prior to observation). Dust extinction for these alternate star formation histories is varied so that the 
UV-continuum slope /3 is equal to the same values (—2.2, —1.5, and —0.8) as for our fiducial SEDs. In general, the f/V-continuum slopes 
/3 we derive from the W colors shown here show little dependence upon the redshift of sources in our dropout samples. Small changes in 
the mean redshift of our selection (i.e., Az ~ 0.1) result in shifts of ~ 0.02, ~ 0.02, ~ 0.01, and ~ 0.05 in the /3's derived for our C/, B, V, 
and i dropout selections (see Appendix A). 
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APPENDIX 

DERIVING THE C/F-CONTINUUM SLOPES FROM THE OBSERVED COLORS 

Here we describe our method for converting the rest-frame W colors we measure for galaxies in our samples to their 
equivalent [/y-continuum slopes. We base these conversions on model SEDs (Bruzual & Chariot 2003) calibrated to 
have ?7T^-continuum slopes /3 of —2.2, —1.5, and —0.8 over the wavelength range '-leOOA to ~2300A (to match the 
wavelength probed by the deep multiwavelength data available for our sample). The model SEDs assume a e~*/'^ star 
formation history, t = 70 Myr, r = 10 Myr, [Z/Z©] = -0.7, and a Salpeter IMF with E{B - V) dust extinction of 
0.0, 0.15, and 0.3 applied using the Calzetti et al. (2000) attenuation law. We have calculated the color of these model 
SEDs as a function of redshift by integrating these SEDs across the filter sensitivity curves. The results are shown in 
Figure [T5l as the solid red lines. We also have calculated the colors expected when the model SED has a perfect power 
law slope /a oc and cuts off at 912A (see black lines on Figure ITC)). 

Then, knowing the effective /? of the model SEDs and the observed colors of these SEDs at the mean redshift of our 
different dropout samples, we can derive a relation between the measured colors and the C/l/-continuum slope. The 
derived relations are 



/3 = 3.04(F - /) - 1.99 



(for U-dropouts) 



(Al) 



30 



13 = 5.30(i - z)- 2.04 (for B-dropouts) 

P = 2.45(z - ( J + H)/2) - 1.98 (for V-dropouts with Jno-band data) 
(3 ^ l.61{z - H) ^ IM (for V-dropouts) 

(3 = 2.47( J -H)- 2.27 (for i-dropouts) 



(A2) 
(A3) 
(A4) 
(A5) 



The mean redshifts assumed here for our U, B, V, and i dropouts in deriving these relations are 2.5, 3.8, 5.0, and 
5.9 (Bouwens et al. 2004; Bouwens et al. 2007). If the mean redshifts of our dropout selections differ from the model 
redshifts by Az ~ 0.1, the derived /3's would change by ^ 0.02, ^ 0.02, ^ 0.01, and ~ 0.05 with respect to those 
quoted here. 

While the conversion formulas above were derived assuming a very plausible star formation history for z ^ 2-6 
galaxies, the actual star formation histories for individual galaxies in our samples likely show considerable variation. 
Some galaxies are likely in the midst of starbursts when we observe them, while other galaxies have likely undergone 
bursts some time in the past and so would have slightly older stellar populations. This results in some degree of 
variation in the shape of the SEDs in our samples - which is above and beyond what we can parametrize using the 
model SEDs above, with a simple one-parameter variation in the overall dust extinction. 

To show the extent to which the observed colors depend on the assumed star formation history, in Figure [T5] we also 
show these colors for a relatively young stellar population (with 30 Myr of constant star formation) and a somewhat 
older stellar population (a model for which star formation proceeded at a constant rate for 100 Myr and then ceased 
for 10 Myr prior to observation). Dust attenuation for these alternate star formation histories is varied so that the 
UV-continuum slope /3 is alternatively —2.2, —1.5, and —0.8. The colors of these model SEDs at the mean redshift of 
our dropout samples are similar to that for our fiducial model (by construction). The only exception is the i — z color 
for our _B-dropout selection which probes a slightly different wavelength range (i.e.. Tabled than the default range 
1650A to 2300A. 

We can attempt to use these model SEDs to estimate the approximate uncertainty in converting from the observed 
colors to J/y-continuum slope /3. Assuming the high-redshift population is composed of an equal mix of the two 
alternate model SEDs and our fiducial model SED and averaging over redshift, we estimate an error A/3 of 0.02, 0.10, 
0.02, and 0.02 in the derived /3's at z 2.5, z ~ 4, z 5, and z ~ 6. 



